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The mathematical problem of the scattering of slow neutrons by chemically bound protons 
is stated in terms of a boundary condition and its equivalent integral equation. By means of 
the latter it is possible to calculate corrections to the results of Fermi which can be considered _ 
as the first approximation to the more accurate equations. The method is applied to free 
protons, and it is shown how the well-known results for this case are obtained from the more 
general equations. This special case allows also an estimate (10~*) of the order of the corrections 
in the case of chemically bound protons in agreement with general estimates in connection 
with Eq. (8.4) of Section 7 of the present paper and the application to harmonic oscillator 
binding made in the succeeding paper jointly with P. R. Zilsel. The corrections discussed 
are greater than the inaccuracies of the boundary condition formulation which are of the 


order 107. 


1. INTRODUCTION 


ERMI! has given a theory of the effect of 

chemical binding of protons on their scat- 
tering cross section for slow neutrons. At the 
time of Fermi’s work there was no necessity of 
great accuracy, and his answer was generally 
considered to be good enough. The advance in 
experimental methods that has taken place since 
makes it desirable to be able to estimate theo- 
retical effects more accurately than previously. 
The whole effect dealt with is of the order of a 
factor 4, and it is natural, therefore, to inquire 
into the accuracy of its estimates. Its knowledge 
is needed in the interpretation of experimental 
material.on neutron-proton scattering in terms 
of effective neutron proton potentials, and the 
comparison of such potentials with corresponding 

* Now at Yale University 


1Enrico Fermi, aay Scient. VII 13 (1936). See 
Eqs. (71)-(103). 


quantities for other pairs of nuclear particles. 
As is well known, such a comparison is essential 
for the field theories of nuclear particles, and it 
has indicated in the past that the specifically 
nuclear interactions are essentially independent 
of the charge. It appeared advisable to be 
reasonably sure of the accuracy of the Fermi 
correction factor in this connection on account 
of its bearing on the fundamental problem of 
forces between nuclear particles. An additional 
incentive has been the large accumulation’ of 
experimental material on the scattering of slow 
neutrons by solids which presents very similar 
mathematical problems. In this case, just as in 
the scattering of neutrons by protons, the nuclear 
properties can be described by quantities a which 
have the significance of intercepts on the axis 
__t Herbert L. Anderson, Enrico Fermi, Leona Marschall, 

, Proc. 7. a Soc. Meeting at Chicago, 
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referring to the neutron-scattering nucleus dis- 
tance in a graph of the product of this distance 
by the radial wave function. The quantity a has 
the significance of the amplitude of the scattered 
- wave per unit intensity of the incident wave in 
a suitable normalization. The interactions of the 
scattered neutron with individual nuclear parti- 
cles, matter in the problem of scattering of slow 
neutrons only insofar as they affect the inter- 
cepts a which are therefore important nuclear 
parameters. The close relationship of the present 
paper to the mathematical problem of deter- 
mining these parameters added to the desirability 
of finding a systematic way for the calculation 
presented below. 

The method consists of two steps: (a) the 
replacement of the wave equation in the region 
of the nuclear potential by a boundary condition 
Eq. (3); (b) the formulation of the boundary 
condition in terms of an integral equation, Eq. 
(5). The physical justification for the boundary 
condition is made in Section 3, and it is found 
that the errors caused by this replacement are 
negligible in comparison with all effects con- 
sidered here. The conversion of the problem to 
an integral equation is made in Section 4. 
Fermi’s approximation is discussed in Sections 5 
and 6, and the reasons for undertaking the 
present calculation are gone into in these sections. 
In Section 7 the expressions for the “‘back- 
ground” are transformed into forms more suit- 
able for applications than the original. In Section 
8, the method is illustrated by considering free 
protons and seeing how the well-known answer 
is verified consistently with the approximations 
made. 

2. NOTATION 


(Xz, Yx, =proton coordinates. 

I, = (x,, Z,) =neutron coordinates. 

U=U(xzs, yx, %)=potential energy of proton 
caused by molecular binding. 

r=r,—r,=(x, y, ) =relative coordinates. 

E,=energy of proton in a molecule in mth 
excited state. 

u,(fz) =characteristic function of proton in state 
E,. 

A, 

A,= 0? /dx,? +0 

A 


’ ¥=wave function of proton and neutron. \ 

g(r) =potential energy of proton in the field of 
the neutron ; assumed to be a function of the 
distance r= |r| only. 

E=energy of system. . 

M=proton mass; the mass of the neutron js 
assumed to be equal to the mass of the proton, 

V=vector operator (gradient) with components 
d/dx, 9/dy, 9/dz. This is distinguished by 
suffixes v, x to indicate replacement of relative 
by neutron or proton coordinates. 

p=distance of the order of 3X10-" cm beyond 
which the function g may be assumed to be 
zero. 

a=intercept of tangent to graph of r times radial 
function for proton-neutron motion plotted 
against an axis of r. The intercept is by 
definition positive if tangent cuts axis of r to 
the left of origin. 

dr = dxdydz =volume element in space x, y, 2. 

dr, =dx,dy,dz,. 

dr, =dx,dy,az,. 

=center of mass coordinates. 

(A)=angular average of any quantity A over 
directions of the vector r. Equal weights to all 
solid angles are given in the average. 


3. THE BOUNDARY CONDITION REPLACING THE 
PROTON-NEUTRON INTERACTION 


The wave equation will be taken to be 
Ey= (h?/2M)(4,+4;) + U+g(r) (1) 


in the notation just explained. It is assumed 
that it is good enough to use one function g(r) 
so that the discussion is applicable only to slow 
neutrons since otherwise one would have to bring 
in phase shifts for p and d waves which may 
correspond in general to other interaction po- 
tentials. The wave function y in the above 
equation is supposed to correspond to a definite 
value of the total spin. Different values of g(r) 
are supposed to be used for the triplet and singlet 
interactions. This means, of course, that the 
representation of an ordinary experimental situ- 
ation involves an analysis into wave functions ¥ 
such as those dealt with in Eq. (1) with different 
values of g(r). Such an analysis is well known 
and is described in most books concerned with 
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the theory of nuclear physics. It will not be gone 
into, therefore, at this stage. 

The eigenfunctions u, describing the stable 
states of a proton in a molecule satisfy the 
equation 


Entin = (— U) Up. (1.1) 


The wave function y can be analyzed into a 
linear combination of the functions u,(r,) with 
coefficients which are functions of the proton- 
neutron distance vector r: 


Y= Lon (Te) (1.2) 


The quantities a, occurring in the above formula 
are taken to be constant while the functions x, 
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absolute position of the neutron with respect to 
the proton. 

It will be shown next by means of a few 
estimates through substitution of this form of ¥ 
into Eq. (1), that for small values of the proton- 
neutron distance r, the difference in the space 
character of different functions x, is negligible. 
This means that in a good approximation one 
will be able to describe the function y as a 
product of a function of r multiplied by a 
function of r,. This approximation will be used 
only for values of r of the order of nuclear 
dimensions. The argument for the separability 
of y into two factors in this range of values of r 
follows. 


depend only on the relative rather than the Substitution of (1.2) into (1) gives 


It is understood here that the x occur only with relative neutron-proton and the u only with proton 
coordinates. For a fixed value of r, integration of the last equation over the proton coordinates 
after multiplication by u,* gives 


1 ‘ 


where the operators p stand for (4/1)¥ and correspond to the usual representation of the momentum. 
It is now clear that if the right side of (2.1) were absent, then inside the range of action of the nuclear 
potential g, one would be dealing with separate equations for each x,, and that there would exist 
solutions in which only one of the a, is distinct from zero. The states u,x, can be considered as 
being uncoupled in this approximation. Further, it is readily seen that if one were to set all the 
E, in these equations equal to zero, the results of calculations on scattering would be unaffected 
provided the left side of Eq. (2.1) is mever used for values of r greater than the range of nuclear forces. 
Within these small distances the effect of the very slight change in the value of the effective kinetic 
energy E—E, —g is small, and the change in the logarithmic derivative of x, caused by E, is negligible 
at r=p. The truth of this point of view may be seen by noting that the effective kinetic energy is 
affected by E, to the extent of, say, 10 ev while its total magnitude is of the order of 10’ ev. The 
effect on the logarithmic derivative rd(rx)/(rx)dr at r=p is* of the order 
oF M fr 
—| - xin. (2.2) 
dEL For WF Jy 
F=rx, (2.3) 


where 5Exin is the change in the effective kinetic energy. This change in the logarithmic derivative 
corresponds to a change in the intercept a in accordance with 


*G. Breit and E. Wigner, Phys. Rev. 49, 519 (1936). See Eq. (29). 
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By means of (2.2) one can transform the right side of the above equation and obtain F: 
ba pta bE 1 
p Exm pF? 
p?/M =Exin. (2.5) 


The factor (1/pF*) So’ F'dr in (2.4) is of the order of } since this is the value of this factor for an 
interaction potential g(r) in the form of a square well. The quantity p*p*/h? is of the order 
(2m)?(A/4)?/A? = (#/2)?=2.5. Here \ is the wave-length of relative proton-neutron motion inside the 
potential well. With these approximate values one obtains by means of (2.4) 


6a (p+a)? 


a pa 


F’dr, ( A) 


where 


For a=20X10-" cm, p=3X10-" cm, this relation becomes 


ba/a = 105E/ Exin, 
and for 5E=10 ev, Exin=10 Mev 
ba/a10-. 


The corresponding change in the scattering cross section is 0.002 percent and is negligible in most 
problems. The magnitudes of the quantities used in the estimates just made have been such as to 
give a too large rather than a too small value for the effect of 5E, and it may be considered justifiable, 
_ therefore, to neglect its effect. 

So far it has been shown that one has a negligible error if one neglects the difference between E 
and E—E, on the left side of Eq. (2.1). It is seen also that the coupling between different a, caused 
by the occurrence of non-vanishing right-hand members in this equation is similarly of little conse- 
quence. These terms are of the same order of magnitude as the terms just discussed and contribute 
approximately the same amount to Ax,. They may usually be neglected, therefore. An exception 


to this general situation occurs if one of the a, is very small compared with some other a, or with | 


at least one of them. If this is the case, the right side of Eq. (2.1) cannot be neglected because it 
may contribute an important part of Ax,. In order that this should take place the ratio a,/a,, has 
to be of the order of 10-*, and this means that one is dealing with an unimportant part of the wave 
function in the expansion (1.2). It is thus seen that for ranges of nuclear forces that are customarily 
assumed it is legitimate to neglect the coupling between the a, as well as the differences between 
the x, provided, of course, one stays in the part of configuration space r <p. 

This situation suggests a reformulation of the problem in terms of a boundary condition at r=0. 
The quantity p is made to approach zero. Simultaneously g(r) is varied in such a way as to keep a 
at its experimental value. As the limit of p=0 is approached, the errors in the cross section assume 
values much smaller than the above estimate of 0.002 percent and finally approach zero. In view 
of the fact that the difference between this idealized and the more realistic original picture is negligible 
for most practical purposes, the limit of very deep and narrow potential wells will be used in what 
follows. This means that as r—0 the limiting form of any x will be taken to be 


(1+a/r)b, 
where b is a constant. Accordingly the limiting form of ¥(r,, rz) for |r,—r.| =r—0 will be taken to be 
limy(r,, rz) =(1+a/r)f(r.). (3) 


In this equation, the function f is a three-dimensional function of the proton coordinates. It is so 
far an unknown function and has to be determined through the requirement that y behave properly 
in other regions than those corresponding to near coincidence of the proton and neutron. 
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4. CONVERSION TO AN INTEGRAL EQUATION 


By means of the boundary condition (3) the problem is simplified through the fact that through 
all of the configuration space with the exception of the three-dimensional region r=0, one may 
require instead of (1) the more manageable equation 


Hw = Ey, (4) 
— (h?/2M)(4,.+4,) + (4.1) 
One has to find solutions of this equation subject to the following requirements : 


where 


(a) Validity of boundary condition (3). 

(b) The correct behavior of y in regions of configuration space corresponding to large values of 
|r,|. In these regions one should require that the wave function be representable as the sum 
of a plane neutron wave with a proton in a definite bound state plus diverging neutron waves 
with the proton in different bound states. 


An explicit construction of such a wave function can be made by means of the quantity 


\ 


K (rs, r,’) u,(r ’ (4.2) 
h? | | 
where 
ke =(2M/h*)(E—E,) (4.3) 


with the understanding that for real x, one makes x,>0, and for pure imaginary x,, ix, is negative. 
This choice of the signs for x, has as a consequence that, in Eq. (4.2), the separate terms contributing 
to K are either diverging waves as is the case for real x«,, or else they are attenuated waves which 
vanish exponentially as: |r,—r,’|—>«. The construction of the kernel K is very similar to that 
employed‘ in the discussion of resonances in nuclear reactions. The especially useful property of K 
is expressed by the equation 


(H)—£) f K(rz, t,; t,’)D(r,’, 1,’)dr,'dr,’ = D(r,, (4.4) 


- where D is a so far arbitrary density function. It will be convenient from now on to employ the 
abbreviation 


K(t,, %')=K(, (4.5) 


‘and to express the wave function y as 


+ f K(, (4.6) 


where Yo(te, r,) is the incident wave 


Here u; is the eigenfunction of the proton in its “‘initial’’ state, and p; is the momentum of the 
neutron in its incident condition. The form (4.6) assures the correet behavior of y for large |r,|. 
The density function D has to be determined next so as to satisfy the requirements at small r. 
Substitution into (1) with the aid of (4.4) shows that 


gv = —D(r,,1,). (4.7) 
One obtains on substitution into (4.6) 


*G. Breit, Phys. Rev. 58, 506 (1940). See Eq. (13.4). 
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This integral equation has been obtained without approximations from the original differentia] 
equation. It does not involve the approximations of the shortness of the range of nuclear forces or 
any of the specializations made in connection with the introduction of the boundary condition for 
small r of the preceding paragraph. One could employ this integral equation if it were necessary to 
take into account a finite range of nuclear forces. Since the arguments of the preceding paragraph 
have shown that the error caused by neglecting the finite range is negligible for ordinary problems, 
the integral equation (4.8) will now be modified so as to eliminate g(r) and to bring in instead the 
intercept a. ' 

The transformations involved in this process are described in Appendix I. It is shown there that 
in the limit of a very short range of nuclear forces 


where f is that of Eq. (3). Substitution into (4.8) gives now 
ixg| 
(1+a/r)f(t2) =Yo(te, Tr) | (r—0). (4.95) 


This equation has been averaged over the directions of the small vector r because otherwise one 
does not obtain a limit independent of the direction in which the point r=0 is approached. The 
reason for averaging over directions is apparent also in the discussion of Eqs. (5.3) to (5.7). The 
_ boundary condition at r-0 is expressed directly by Eq. (4.95). It is inconvenient to leave the con- 
dition in this form, however, because for r=0, the left and the right sides become infinite. Substi- 
tution for f(r,)/|r.—r,| of its expansion in terms of the u,(r,) gives a neater result, viz. 


with 
Bry, 1) u(te) | f(r’) dr’. (5’) 


In this form every term approaches a finite limit as r= |r,—r,|—>0. Equations (4.95), (5) amount 


to three-dimensional integral equations on the function f(r,), and the problem is now reduced to 


the solution of an integral equation. 

It follows from the way in which Eq. (5) has been derived that the coefficient of a@ on its left 
side should be finite. Otherwise the equation could not hold because all other terms in it are obviously 
finite. This fact may be verified as follows. One finds by direct calculation 


eine! r’—ry| h 


2 
dr’ = 5.1 


It will now be noted that the multiplication by U which is contained in Hy—E on >, on the left 
side of Eq. (5.1) cannot introduce a 6 function unless the 6 function should be a part of the Dy. 
It is clear, however, that for a regular f the 6 function is not a factor in f. The 6 function on the 
right side of Eq. (5.1) has to arise therefore through the application of the Laplacian operators in Hp. 
The combination of these that occurs is . 

h? 

——(A4,+4,) = ——A-— 
2M M 


(Ho—E) Xs tlt) f f(t") 


A. (5.2) 


where the last term represents the Laplacian with respect to the center of mass coordinates of the 
proton and neutron. Only the first term on the right side of Eq. (5.2) has anything to do with the 
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right side of Eq. (5.1) so that the coefficient of 1/r in 2°, can be readily ascertained. It is thus found 


on Dd. of (5.1) =1/2r+---, 

where the remaining terms indicated by dots do not give rise to the singularity. Comparison with 
Eqs. (4.95), (5), (5.1) shows now that the coefficient 2 in the numerator under the integral sign in 
Eq. (5’) is just right to cancel the singularity in 1/r. That this is the case will also be verified in 
the course of applications of the integral equation. 

The integral equation (5) will now be derived without the use of the contraction of range pro- 
cedure. The employment of the interaction potential g(r) will be replaced by the requirement that 
the s part of the wave function is of the form (1+a/r)f everywhere except at r =0. The wave function 
y can be analyzed in spherical harmonics Y, of the angular coordinates of r as follows: 


Ure, 1) + eu(R, 1) (5.3) 


Here r, R are the relative coordinates and the coordinates of the center of mass, the functions gz 
are regular everywhere, while © is not regular at r=0. There is no assumption involved in expressing 
y in the form (5.3). For any fixed R, y is a function of r only, and if r is fixed in addition, then 
one has a function of the angular coordinates of r which can be expanded in terms of the spherical 
harmonics Y,. Averaging Eq. (5.3) over angles one finds, of course, 


(W(R—1/2, R+1/2))=4(R, 7). (5.31) 


It has been assumed that the functions g, (L=1, 2, ---) are regular everywhere so that as r is 
made to go to zero, the limit of any gz is the same independently of the direction from which the 
point r=0 is approached. The function #9, on the other hand, is assumed to be such that 


lim[ 40(R, r) — &(R)(1+a/r)]=0. (5.32) 


The assumptions just made about py and ¢ are what replaces the employment of g(r) in the original 
wave equation (1). In Eq. (5.32) the function #(R) is independent of r as is indicated by the notation 
used. 
It must now be explained that to a sufficient approximation for small values of r, 


(W(R—r1/2, = (te, te (5.33) 
In fact, according to Eq. (5.3), | 


Since there is only s scattering, the functions ¢g, for small r vanish as r¥ or faster. The occurrence 
of /2 in the combination r,—r/2 in g, brings in a dependence on the direction of r in addition 
to that included in Y,. Expanding the gz in Taylor’s series, one obtains terms vanishing at least 
as fast as r*. Such terms will not matter in the statement of a boundary condition such as (5.32) 
because they disappear from y as r approaches zero. One is concerned, therefore, only with 


($0(t,—1/2, r)) = 
= (1+4a/r)((1,) + (x*/8) (0°@/dx,") + +) (5.35) 
=(1+4a/r)[®(t,) +(7°/24) +--+]. 


In the condition expressed by Eq. (5.32), it does not matter, therefore, whether one averages over 
angles with fixed center of mass coordinates or whether one keeps the neutron coordinates fixed, 
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averages over angles of r, and then makes r, approach R. The same holds for r,. Combining this 
result with the consideration of the y, which was mentioned right after Eq. (5.34), one may restate 
Eq. (5.32) in the form 


lim[ (1+a/r) (5.4) 
_ If now y is defined by 
Tr) =Yolts, Tr) (5.5) 
where >", is as in Eq. (4.9), then simply as a matter of rearranging terms 
(5.51) 


with B given by Eq. (S’). Averaging (5.51) over angles in order to substitute into Eq. (5.4) one 
has to consider first of all the term Wo. One has for it 


so that in (5.4) one may replace this term by yo(r,, rf). Substituting all terms that matter into 


(5.4) one has finally 
(1+a/r)®(r,) tx) (5.6) 


where O(r) designates terms of order r. Identifying coefficients of equal powers of r in Eq. (5.6) 
one has 

&=f; =Yo(tr, tx) +a(B),—o. (5.7) 
These two equations are equivalent to Eq. (5). 

It has thus been shown that the boundary condition statements of Eqs. (5.32), (5.4) require the 
validity of the integral equation (5) if it is assumed that y can be expressed in the form obtained 
by substituting for the integral on the right side of Eq. (4.8) its limiting form as in Eq. (4.9). The 
limiting contraction process has dropped out of the picture. Its only remnant is the suggestion of 
a possible form for ¥, and this suggestion is not necessary since one can show that (B),.o is a well 
defined quantity as will be seen in connection with Eq. (6.1). 


5. COMPARISON WITH FERMI’S EQUATION 


In the present paragraph a comparison with Fermi’s equation! will be made from the point of 
view of expressing the quantities occurring in one in terms of those occurring in the other. The 
numerical comparison is postponed until formulas for the evaluation of quantities occurring in 
Eq. (5.1) are worked out. 

Combining Eqs. (4.8), (4.9) one obtains 


|r.’—r,| 
It follows from this and Eq. (5.1) that 
(Ho—E)y = (6.1) 
This relation has some similarity to Fermi’s equation : | ; 
= (ah? (6.2) 


where 7 is the average of y over positions of the neutron within a sphere of radius R for fixed proton 
positions. The radius R is supposed to be large in comparison with nuclear dimensions and small in 
comparison with the size of the molecule. The delta function on the right side of Eq. (6.2) is first 
introduced by Fermi in such a way as to have a constant value within the sphere of radius R and 
so that its integral through the spherical volume is unity. In the applications of the equation Fermi 
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the function 5g by Dirac’s 6 function, and the function ¥ on the right side of his equation 
he replaces by the incident wave in the spirit of Born’s approximation. 
Comparison of the exact equation (6.1) with Fermi’s approximation (6.2) shows some points of 
similarity and also an obvious difference. Eq. (3) shows that f(r,) is approximately Fermi’s : ) 


f(t.) (6.3) 


it follows, therefore, that if one is interested only in the scattering, then, to within the limits of 
validity of (6.3), the Eqs. (6.1) and (6.2) are nearly equivalent because for the calculation of the 
scattering the distinction between y and 7 on the left side of the two equations is immaterial. One 
may look, therefore, at Eq. (6.1) as a substantiation of the approximate validity of Fermi’s method 
of calculation. One sees further that if Fermi’s equation is used not exactly but in Born’s first ap- 
ximation and with the replacement of 5g by Dirac’s 6, then the only error involved is that of 
approximating f(r,) by the incident wave. Equation (5) shows that f(r,) is only approximately 
equal to Yo and that more accurately one should add to yo the term in a on the left side of (5). This 
term will be referred to as the “background” for the following reason. It arose as the difference 
between the term in 2a on the right side of Eq. (4.95) and the quantity af/r. The term in 2a is, 
according to Eqs. (4.8), (4.9), equal to ~—y» and represents, therefore, the effect of the diverging 
waves emanating from the source density f in the vicinity of r=0. A part of this term represents 
af/r. The excess over this is the “background” term. 
Disregarding the “‘background” term there is no difference between the results obtained here 
and Fermi’s method in Born’s first approximation. The validity of Fermi’s result is, on the other 
hand, no better than that of Born’s first approximation as may be seen by noting that 


1 +") = ( 1 (6.4) 


and that calculation with Fermi’s Eq. (6.2) in the second Born's approximation would bring in 
terms having an order of magnitude relative to the zeroth of a divided by the linear molecular 
dimensions. 
6. SIMPLIFIED EXPLANATION OF FERMI’S EQUATION 
The end result of Fermi’s averaging of the wave equation is the replacement of the nuclear inter- 
action potential g(r) by an interaction energy 


ah? 2 


Vr = =— 
aM R 


(r<R). 


Vr=0 (r>R) 


For R slightly greater than p this quantity is large in the region r<R. The introduction of Vr 
is then not of much use and would even lead to erroneous results for scattering if one used Fermi’s 
equation (6.2) exactly. In fact, the law in accordance with which Vr and R should be varied together 
in order to obtain equivalence for scattering is much more closely VrR?=constant than V,R® 
=constant, and it is the latter way of varying Vr with R that follows from Eq. (7). This way is 
therefore not right, and one should not apply Fermi’s equation in this manner. It is thus seen also 
that Fermi’s equation cannot give even approximately correct results for scattering if one replaces 
the function 5g by Dirac’s 6 function and calculates with the resultant equation exactly. The latter 
replacement corresponds to approaching the limit of R=0 in Eq. (7) which corresponds to varying 


the internal phase as 
M 3\a!| 
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which gives for F of Eq. (2.3) a limiting behavior for R<a 


which is indeterminate for R-»0 while the limiting value of this quantity should be 


(—) - R R ( ) 

~—{1-— }. 
(= L(r+a)dr “neo a ‘a 
These two results are obviously inconsistent, indicating that one should only use Dirac’s delta 
function in Fermi’s equation for purposes of calculating the scattering to the first order of Born’s 
method, but that an exact solution of this equation would lead to incorrect results. For this reason 
the explanation in the present section of the paper is concerned not with the equation which Ferm 
recommends for solution by the Born procedure but with the form (6.2) in which the dz function 
is employed, A simplified presentation of Fermi’s equation is taken up next. 

The radius R can be made several times a, and the potential energy Vr can thereby be made 
small compared with the nuclear interaction potential g(r) and yet large in comparison with the 
molecular potential energy acting on the proton. By making R sufficiently large one can make 
the modification of the wave function due to the action of Vr quite gentle and keep the term of 
ay/R small compared with y. Born’s first approximation can be thus made to be a good approxi- 
mation for Eq. (6.2). The margin of available orders of magnitude is not exceedingly great for this 


_ choice of R but is nevertheless sufficient to make this point of view have good sense. Thus, e.g., 


one can make 
R=20(4/(Mm) ic) =1.8X10-" cm, 


which is reasonably large compared with the maximum value of a that comes under consideration, 
vis. a=2X10-" cm. The averaging sphere is still small in comparison with the molecular vibration 
amplitude which is of the order 10-® cm. This choice gives 


3a me* 1 510000 


r> 


R 20° 3 400 


The kinetic and potential energies of molecular motion are seen to be small compared with V, 
and can be neglected in the region r<R. Under these conditions the dependence of y on r inside 
r<R is expressible as the product of a nearly constant factor and the function x. The dependence 
of the latter on r is for an attractive Vr, i.e., a>0 


F=rx=const. Xsin[(M| Vr| /h*)#r] 


(dF_/F M (7.1) 
rd F_ ~1—-| — =1-—, A 
R MRK R 

which is precisely the value needed to join on smoothly to the function F=r+a for r>R. In fact 
for the latter 


ev = —433 ev. 


so that 


(rd F,./ nar =R/(R+a)=1—a/R+:::. (7.2) 


If Vr is repulsive rather than attractive, the sign of the terms in 1/R in both of the above expres- 
sions is changed, and the internal function can be still joined smoothly to the external one. It is 
seen, thus, that the potential energy Vr is such that for Ra the solution of Eq. (6. 2) gives the 
correct value of rdx/xdr. Within the limitations already discussed this equation must give, cerca 


the solution of the problem. 
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It has been noted above that the margin available for the choice of R is not very great. The 
question naturally arises as to the order of precision of calculations that assume the validity of 
Eq. (6.2) which has been seen to imply the possibility of choosing R so that it is both large in com- 
parison with the range of nuclear forces and small as compared with molecular dimensions. 


7. TRANSFORMATIONS OF THE BACKGROUND 


In order to estimate the accuracy of the calculation, one has to know the coefficient of @ in Eq. (5). 
For this purpose it will be transformed into a different form. By means of Eq. (1:1) one obtains 


and one further has for any regular g(r) 


eur 
f —4n9(0). (8.1) 
o(t) =f(r)u.*(r) 
and combining Eqs. (8), (8.1) one obtains 
Ly f = (2rh?/M)f(r.)>. terms. (8.2) 


By means of a similar transformation one finds 


Dd. f Us = terms. (8.3) 
r,/—f, 


In Eqs. (8.2), (8.3) the first term is the same. This term approaches f(t.)/r for small r. If, therefore, 
one subtracts (8.3) from (8.2), the result is free of the term in 1/r. It is thus found that 


(ts, Tx) E, E, J |r—r,| 


h? 
x f(t)u.*(r) U(r) f(r): Vu.(r)+u,* f(r) ] far. (8.4) 


The behavior of the first term in this formula will now be discussed. One finds by straightforward 
calculation that 


[4,—(2M/h*) U(r.) |X = 
X 
It follows, therefore, that for small r = |r, —r,| X behaves as 1/r. For large r and fixed r, the quantity 


X vanishes much more rapidly than 1/r on account of the rapid decrease of the ,(r,). In + 4) 
the first term contains on performing the summation contributions of the order 


(M/2xh*) Ef/|r—r,|? 


where 


_ for small value of |r—r,| modified by the presence of the exponentials. This term has, therefore, 


tion 
fact 
1.2) 
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the approximate magnitude \ 
2M R E f (tx) 


(8.5) 


where & is a mean length at which the changes in f and in the exponentials containing x, decrease 
the integrand appreciably in comparison with the approximation just mentioned. This length js 
roughly of the order of 1/x, where x is a mean value of x, and the length is, therefore, of molecular 
dimensions. The energy of the neutron enters, therefore, in comparison with the centrifugal barrier 
for molecular dimensions, and besides one, has the factor {/R. Therefore, one expects the order of 
magnitude of the effect of the background on the cross section to be 10-* cm/10-* cm=0.001, 
which amounts to a tenth of a percent. 
Another convenient form is found by introducing the quantity 


where Ep is the lowest E,. One verifies readily that 

U(re)]} F=9F/dr, (8.7) 
and one obtains 

2M 
s=0 

where 


which is useful for the evaluation of B (Cf. Appendix II). This transformation is especially con- 
venient for E= Ep. It can be employed also for other values of E..In applying Eq. (8.8), one repre- 
sents the collection of terms due to the —1 in the numerator of Eq. (5’) as 

1 1 exp(—?7*/47) 


dr. 8. 


It is found then that in the limit of r=0, the two definite integrals over + entering the representation 
of B combine to give a convergent integral. In evaluating B for E> Ep» the number a in (8.8) must 
be chosen so as to make E, <E<E,,;. The terms >> ,.0’f, represent the effect of the diverging waves. 


8 FREE PROTONS 


The equations will be illustrated first in the case of free protons. There might be reasonable 
doubt regarding their applicability because the free protons recoil and the proton waves have 
sources of diverging waves on a par with the neutron waves. These sources are not any present 
- in Eq. (6), however. 

The protons will be quantized in a aeanemented volume with a periodic hisueiie condition. 
The wave function yo will be taken to be 


corresponding to the collision of neutron and proton waves of unit amplitude. The momentum 
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of the incident neutrons is po. The proper functions for the protons are 
u,(tx)= k,=p,/h, (9.1) 


where V is the fundamental volume. In the usual manner the summation over s can be replaced by 
V f dk, /2x*. (9.2) 

One has thus for the sum that is multiplied by 2a on the right of Eq. (6) 
fiv-si- exp ker’ + (bt | (9.21) 


where in accordance with the convention that has been followed for « one must set 


| (ko >k) (9.22) 


The quantity needed for the integral equation is 


1 sinkr 


On introducing =r’ —r, and integrating over directions of r’’ one finds 


f sin(kr) sin(|ko—k|r’’) 
0 


and integration over directions of k gives 


Qetkory 2 
f f (1/r’”) sin(kr) sin (Rer”) sin(ker”) exp (9.26) 


It will be noted that the integrals over r’’ converge and that the integrations can be performed 
first over by straightforward calculation. 


LP. f /r'’) sinkor” sinkr’’ exp[i(ko? — k*)* 
0 


= (i/4) f (1/r) (9.27) 
0 


with 
Ai=V+ho—k; As=vV —kotk; 
As=V+khot+k; As=v —kho—k; (9.28) 
and 
Vv = (ko? —k?*)!. 


Contributions to the imaginary part are present only if ko>k. For this reason Ai, As, As are positive, 
while A, is negative in Eq. (9.27). The right side of (9.27) becomes wi/4 and 


1 —coskor 
IPAS.) = (9.3) 
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Similarly one obtains | \ 


R.P. i kr | = 


and in the limit of r=0, and r,=r, one has from Eggs. (9.3), (9.4) 
| a(B) pmo = takoe™** /2, (9.5) 


Here use was made of the fact that one may keep either r, or r, fixed in the averaging over directions 
of r as has been explained in connection with Eq. (5.34). One may replace, therefore, the r, in Eqs. 
(9.3), (9.4) by re. 

In the above hitless one has introduced into the calculation of B the approximate value of 
f&vWo which follows from (5) if one neglects a(B), i.e., if one solves (5) by iteration. On account of 
Eq. (9.5) the corrected value of f is 


f(t (9.6) 
The asymptotic behavior of the angular average over directions of r is 
(9.7) 
which can be ennipesed with the exact value following from 
according to which = 
sin (koa /2) sinkoa/2 


1+ (1 /r+iky/2) = +coskea/ +ikoa/2)(1+a/r). (9.8’) 
ho/2 r 


koa/2 


. It is seen from Eq. (9.8’) that the wave function is not quite of the form 1+a/r. The deviations 
from this form are seen, however, to be of the order k,’a?, which is a very insignificant factor for 
thermal neutrons. 

It is satisfying that to the order koa there is agreement between (9. 8’) and (9.7). While even the 
terms in ikoa/2 affect the cross section only by the insignificant factor 1+:’a*/4, it is satisfactory 
to have the phase of the scattered wave correct to the order koa/2. 

The agreement between (9.7) and (9.8) is seen to be no worse than the assumption that the factor 
1+<a/r represents correctly the ‘variation with r of the s part of y at small r. The exact dependence 
is of the type 

tanka a 
(inka coskr +coska sinkr) ~costal 1 =| 
r 


The deviations from the assumed form are themselves representable by correction factors 1+k’a/3, 
and there is no apparent reason for expecting the agreement to be closer than that found. 

If one makes a calculation for the cross section along standard lines with yo in place of f in Eq. (6), 
then one finds without difficulty that the differential cross section per unit solid angle in any given 
direction when corrected for the factor 1+ ’a?/4 is 


w= (40° + ko?) /po. 
The last correction goes outside of the limits of accuracy of the method and does not represent an 
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improvement over the simpler result 
4a*p,/ Po, 


since this correction is given by the factor 
sin*koa/2. 


If all the corrections obtainable from the equations were as small as that just discussed, there would 
be no purpose in presenting the equations. It is found, however, that the actual correction factor 
involves the ratio of a to molecular dimensions to the eet rather than to the second power and is 
not negligible, therefore. 

A general idea of what one may expect can be obtained by considering the background for free 
protons as made of two parts, one of which comes from the region k, >ko, and the other from k, <p. 
The evaluation can be carried out by means of Eq. (9.4) by replacing the integral by 


k+ko _—coskor sinkor 
f sinkr log — 2h) (9.8 
ko 


with 
C=0.5772---= Euler-Mascheroni constant 


Si(x) = f Cils) = f du. 


By means of the above Eq. (9.9), one finds on substitution into Eq. (9.4) and on passing to the 
limit of 


(10) 
k>ko 
and one has similarly : 
RPK > (10.1) 
k<ko 4 


It is seen that the singularity in 1/r is contributed by the region of high k as could be expected since 
the representation of the 6 function in Eq. (6.1) is impossible without the extension of the sum over 
the eigenfunctions to those of infinitesimal wave-length. It is seen also that there are two real contri- 
butions to a(B), one of amount 


—ak)/x=part of (B),-o for k> ko, (10.2) 


and another of equal and opposite amount which comes from k<ko. The effect of each of these 
parts on f is of the order of a fractional change of amount 2a/Ao, where Xo is the wave-length of the 
incident neutron which is of the order 10-* for 1 ev neutrons, and 10-* for 100 ev neutrons. When 
the protons are chemically bound, the balance of the contributions from k>ko and k<kp is upset, 
and there is a residual effect which can be expected to give changes in the scattering cross section 
of the order of tenths of a percent, since the molecular energies are of the order of 0.5 ev or 0.2 ev. 

The evaluation of the background for neutrons having very small energies incident on harmonically 
bound protons is made in the next paper in this issue.* The results agree with the estimates made 
here in Section 7 above. 


*G. Breit and P. R. Zilsel, Phys. Rev. 71, 232 (1947). 
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APPENDIX I 
Derivation of Eq. (4.9) 
Substitution of K into the left side of Eq. (4.9) by means of Eqs. (4.2), (4.5) gives 


f K( ’ "\dr,'dr,’ (Fe dr,’ (I, 1) 


where the replacement 
=f(te)x(r) (I, 2) 


has also been made. The latter involves the harmless approximations which have been justified in 
connection with Eqs. (1.2), (2.1), and (3). This discussion showed that to a very good approximation 
[ ) =0 | 
r) Ix(r) =0. 
re g(r) Ix (I, 3) 


As the range of nuclear forces p contracts in the process discussed for Eq. (3), oné can make it small 
enough, for any pre-assigned r,’, r,, so as to have 


|r,’—r1,|>p. (I, 4) 
When this condition is satisfied, one can set 
etteltr—te’| 
where 
I= (1, 4.) 
and the employment of Eq. (I, 3) in the above formula gives 
kta 
(I, 4.3) 


In the calculation just described a appears as the value of —r*dx/dr at the outer boundary of the 
region where g(r) #0. Substitution of (I, 4.3) into (I, 4.2), of the latter into (I, 4.1), and of the 
result into (I, 1) yields Eq. (4.9) of the text. 

In the above discussion no account was taken of the region in which (I, 4) is not fulfilled. As p 
approaches zero, the volume of the part of the space r,’ in which (I, 4) does not hold also approaches 
zero, and the contribution to the sum on the right side of Eq. (4.9) due to an individual term is 
of the order 

(I, 4.4) 


This is true only for terms with sufficiently low s. For high s, the strongly negative numbers ix, 
make the contributions smaller. The series consists of terms like (I, 4.4) up to the point where |«,| 
is of the order 1/p. The sum on the right side of Eq. (4.9) contains thus 


lxs|~1/p 


| 

| 

| 
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which is an approximation to the 6 function and differs from the 6 function by being distributed 
through a volume of approximate linear dimensions p. The whole sum on the right of Eq. (4.9) is 
thus of the order f(r-)/p, if r=|r,—r,| is of the order p on account of (I, 4.4), (I, 4.5), which has 
to be true in order to cancel the term in 1/r on the left side of Eq. (4.95). In the description of the 
behaviour of the sum which has just been considered, the quantity p has no direct significance as 
the range of nuclear force, but can be made equal to it if desired. It is thus seen that the contribution 
to the sum owing to the region |r,’—r,| is quite large for |r,—r,| <p and cannot be said to be 
negligible everywhere. It is also seen, however, that on account of (1, 4.5), this contribution is very 
much smaller if |r,—r,| >p, because the D function in (I, 4.5) is concentrated mainly in a region 
having linear dimensions p. . 

In spite of the fact that the contributions to the sum in (4.9), caused by the region inside a sphere 
of radius p, are responsible for the term f(r,)/p, when |r,—r,| =p, the contraction procedure makes 
it possible to make p small enough so that for a given |r,—r,| one can satisfy (I, 4). 


APPENDIX II 
Derivation of Eq. (8.8) 
For real values of ix,, one can represent the quantity f, which is defined by Eq..(8.9) as 


fam f exp { |r,—1’|?/47— | ke|2r}dr. (II, 1) 


This formula can be verified as follows. The factor in (8.9) which multiplies the product of u, and 
u,* will be temporarily denoted by g,. One can define g, by 


(4,+«,*)g.=0, (II, 1.1) 
and the requirements: (a) g.—0 for r, ; (b) g, approaches |r,—r’|—! as a limit for |r,—1’| 30. 
It is seen that ' 


4 exp[— |?/47— | (II, 1.2) 
0 


has the property (II, 1.1) of g, 
o ar 


and is, therefore, zero unless |r,—r’| =0. In the latter case the right side of (II, 1.3) is infinite. 
It is also clear from the form of (II, 1.2) that hk, has the same asymptotic forms at small and at 
large distances as g,. It follows that 

g.=h, (II, 1.4) 


as can be verified by integration of (II, 1.2). Substituting (II, 1.2) for g, in Eq. (8.9) of the text 
one obtains Eq. (II, 1) for f,. This transformation applies only for x,2<0 because if x, is real, then f, 
is an outgoing wave. Equation (8.8) now follows by straightforward substitution of (II, 1) for 
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\ 


The cross section for scattering of very slow neutrons by elastically bound protons is calcu- 
lated by means of the method described in the preceding paper. The results are expressed 
in terms ef a correction factor to be applied to Fermi’s factor 4. According to Eq. (5) of the 
present paper this correction factor is 1+16aJ/x+ with J defined by Eq. (6) which is suitable 
for numerical work. The results of numerical calculations for the effects of anisotropy in 
molecules with axial symmetry are presented in Fig. 1 which has to be used with Eqs. (6.2), 
(6.3) of the text. For scattering by the H:O molecule the correction to be applied to Fermi’s 
result is only ~.3 percent of the final value. Fermi’s method is very accurate in this case. 


1. INTRODUCTION AND NOTATION 


HE method of the preceding paper’ is here applied to the calculation of the scattering cross 
section for neutrons having zero energy from protons bound in molecules by simple harmonic 
oscillator forces. The F transformation of the background term! is used to reduce the computation 
of the background to the evaluation of a definite integral by numerical quadrature. It is found, in 
agreement with order of magnitude considerations discussed in the preceding paper, that Fermi’s 
factor 4, representing the increase in the scattering cross section caused by chemical binding, is 
accurate to within a few tenths of a percent, and the value of the correction factor is worked out 
in the last section of this paper. The notation is as in' with the addition of the following symbols: 


a=constant having dimensions of length defined so that radial wave function is of the form 
const.(1++-a/r) for small proton-neutron distances r. 
@z, ay, a,=elastic binding constants of proton in molecule. 
n=y/ay; 
H,.(&) = (—)*e(d/dt)"e-® = Hermite polynomial of degree n. 
; gy=exp(—27/a,*) ; 
etc. 
The potential energy binding the proton in the molecule is taken as 
h? Jur 
U(r.) =—{ —+—+— }, (1) 
2M 


as! 


for which the characteristic frequencies are 


2nv.=h/Ma,’, etc., (1.1) 
so that the energy levels are given by 
: (2my+1)/ay?+ (2m. +1) /2M (1.2) 
and the proton eigenfunctions by? 
u(r) = exp(—x*/2a,*) (1.3) 


* Now at Yale University. 

1G. Breit, Phys Rev 71, 215 (1947). This will be referred to as I in the text. Formulas occurring in I will be referred 
to by prefixing the roman numeral I ahead of the number of the formula. 

* In this paper the symbol II always stands for the product over ~, ¥, 2. 
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with 
H,() =Hermite polynomial. (1.4) 


The following notation will be used 
E=x/az, n=y/ay, (1.5) 
The complete wave function of the system is given by Eq. (I, 6) with f given by (I, 5) and B by 
(I, 5’). 
2. CALCULATION OF THE BACKGROUND 
Equation (1, 5’) can be evaluated by means of the transformation (I, 8.8). The firstapproximation 
to f is Yo, and the solution of Eq. (I, 5) by iteration leads, therefore, to an integral in which f is 
replaced by yo in (I, 5’). In the present case, with E= Ey and up given by (1.3), one can replace ¥ by 


(2) 
Putting H,,(&) =(—)" exp(é)(d"/dé") exp(—) in (1.3) and substituting into Eq. (I, 8.6) one has 
+#9)/2 
F(r, r’) = Ding, ny, ms “(0 (2.1) 


The summations in (2.1) can be performed by expressing the exponentials containing ——£? by 
means of Fourier-integrals. One obtains 


1 
exp(— f exp(iwt —w*/4)dw, 
and thus 
1 @ 
exp(— #— #7) =— f f expLi(we-+w! /A dade’. 
4n J_.. . 


Puttin 
gz=exp(—2r/a,*), gy=exp(—2r/a,?), (2.2) 


one obtains the auxiliary formula 


exp(— &2)/2"(n!) = Fe 


1 2 2 2 n 
-—f (g/2) (—ww’)" exp[i(wt+w’ — (2.3) 


=n! 
1 


The subscript x has been dropped for a, m, and g so as to simplify the notation. 
The integrations over w, w’ can be carried out by transforming the homogeneous quadratic form 
in w, w’ in the exponent of the last form of (2.3) to a diagonal form. It is thus found that 


F,=(1—g.*)+ exp{ &?— ]/(1—g,") }. (2.4) 


Comparing (2.1) and (2.3) one sees that 
F=]] exp[(#+ 
exp {(#+ &)/2—(#+ } 


so that 
(3) 


F(r, 


| 
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It follows from Eq. (I, 8.8), on noting that E= Ep, and thus setting o=0, that \ 
2 1 

f exp { } _ ote 
0 2 
exp(— 
x— | (3.1) 


where 1/|r,—r1’| is represented as 
(1/2) f Jar. 
0 


Putting (3.1) and (2) into (I, 5’) and performing the integration over r’ it is found that 


* B(r,, = (8/4) f(rs) f {II (47.)+ exp { — —exp(— (3.2) 
0 


with 
(3.25) 


The limit of (3.2) as r,r, is independent of the angle (r,, r.) so that the averaging process required 
for Eq. (I, 5) reduces to 


B(te, = (8/x)f(ts) f exp { x} — dr. (4) 
0 


In the calculation leading to this equation it was assumed that f'can be replaced by a constant 
multiple of uo. This is justifiable only if the molecule is originally in its lowest state, if the neutrons 
have a very low energy and if a is sufficiently small to make aB small compared with f. The latter 
assumption is justified a posteriori by the smallenss of the effect of B. 


3. THE SCATTERING CROSS SECTION 
For zero energy neutrons the amplitude of the scattered wave is proportional to 


f 


f is given by (I, 5) with B from (4), in which f is taken as f=yo= wo. 
Thus, 


= f fi Jar (80/2!) ay a2) (5) 


The square of this factor, i.e., the quantity 1+16aJ/r', is the correction factor for the cross section. 


The quantity J occurring in this formula is given by 
0 


as is found by performing the integration over r’. 
The integral occurring in Eq. (6) has been evaluated numerically in a few cases,‘and the results 
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Fic. 1. Plot of atl (a) US. Curve Ais drawn by MEAS 
through points computed by means of Eqs. (6), 
(6.1). Curves B are the asymptotic forms given Nite P /) 
by Eqs. (6.3), (6.4). The ordinate of the graph \. 
factor of Eq. y means s. (6.1), (6. 
with @ defined by Eq. (6.3). 


\ 


4 


are plotted in Fig. 1. In this figure the ordinate is the quantity 
a'I(a) =atJ(1, 1, a). (6.1) 


The asymptotic forms of this quantity for large and for small @ are indicated on the graph, and 
their derivation is described in the Appendix. By means of the graph one can evaluate J for any 


harmonic oscillator with axial symmetry by means of the following formula 
J (az, &2, &s) = (6.2) 
The ordinate plotted is 
(a) (az, a2, a2); a=a,/ae. (6.3) 


The asymptotic forms of J are 


J (az, Oz, Os) (a>1), 
where 
h=(1/16) f (04-1) (6.5) 
(1/16) f (07-1) (6.6) 
(6.2) 


The two asymptotic branches extrapolate to nearly the same value for an harmonic oscillator with 
spherical symmetry and agree rather well with values computed directly by means of Eqs. (6), (6.1). 


4. APPLICATION TO THE WATER MOLECULE 


In order to apply the present theory, the mass of the proton will be considered to be small in 
comparison with the mass of the molecule. One has then 


6.39 X10-" cm 
[(hv)ev 


where the subscript ev indicates that the electron volt is the unit of energy. For H,O Slater? lists 
the following frequencies: 


hy (ev)= 0.443 0.462 0.248 
*J. C. Slater, Introduction to Chemical Physics (McGraw-Hill Book Company, Inc., New York, 1939), p. 146. 
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which correspond to \ 
a@ (cm)=9.6X10-" 9.4X10-" 1.28 10-*. 


The values given by Herzberg‘ are in essential agreement with the above. The lower frequency and, 
therefore, the higher a corresponds to a vibration roughly at right angles to the other two. The value 
of a,/a is 1.35. According to Fig. 1 the anisotropy has no important effect in this case. The values 
of a for singlet and triplet scattering are a;=1.92X10-" cm; ag= —0.584X10-" cm. For singlet 
scattering the correction factor is 


1.92 10-2 
1+0.82-—_——= 1.0015. 
1.06 10-* 


The square of this factor matters for the scattered intensity. There is thus a correction of about 
0.3 percent for the expected scattered intensity in the singlet state. Since slow neutrons are mostly 
scattered in the singlet state, this is the approximate magnitude of the correction for very slow 
neutrons. 
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APPENDIX 
Asymptotic Forms for I(a) 
The consideration will be carried out for a1. It follows from Eqs. (6), (6.1) that 


The integral will be represented as 
with 


I(a) =A (a) +B(a), (A, 1.1) 


0 


where the integrand is the same as in (A, 1). The number N is chosen so as to satisfy the requirement 
N>1, No®<1. (A, 1.3) 
In the limit a=0 both of these requirements can be satisfied perfectly in the sense that N can be 


made to approach while Na* approaches zero. Integral A(a) can be transformed through the 


introduction of the variable ' 
§=2r/a. (A, 1.4) 


On account of the second condition (A, 1.3) the quantity e~** can be expanded by Taylor's series, 
and the contributions te the integral form a rapidly converging series. Keeping the first two terms 
of the expansion one has 


A(a) = (1/16) (1 —e-*) 4-1} (A, 2) 
0 


The discussion of the contribution of the term in a’£/4 will be postponed. It will be shown later 
that its contribution to A approaches zero as a approaches zero. Only terms that do not approach 


‘G. Herzberg, Infrared and Raman Spectra (D. Van Nostrand Company, Inc., New York, 1945), p. 161. 
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zero with a will be kept. One has thus 


The first integral will be kept in the above form. The second can be estimated as follows 
1 
-— 


(A, 2.2) 


8a(2N)* 


The right side of (A, 2.2) becomes infinite for a=0 on account of the second condition (A, 1.3). 
The omitted terms, however, are of order N-4a~!. They can be made to approach zero by requiring 


N*o?>1. (A, 2.3) 


This means that N should be made to approach © more strongly than a which can be done con- 
sistently with the second condition (A, 1.3). Any dependence between a~* and a+ is satisfactory. 
With this additional restriction on N it suffices to consider only the right side of (A, 2.2). 

In the evaluation of B one makes use of the smallness of a?/r<1/N. At the lower limit of integration 


[-+=a- (-2 (A, 3) 


The omission of further terms here is justifiable because B turns out to be of the order 1/(Na*)* so 
that (A, 2.3) will make negligible the term in 1/(4N) which has just been omitted. One obtains by 


means of (A, 3) 


The right side of the above relation contains J;"**. Its contribution is estimated by expanding the 
brace in a Taylor series. It is found to be —(Na*)*/16 which vanishes on account of the second con- 
dition (A, 1.3). 

Adding B to A the last term in (A, 3.1) is canceled by (A, 2.2). The integral J, of Eq. (6.5) of the 
text arises from (A, 2.1), the integral J, of Eq. (6.6) from (A, 3.1) and one obtains 


I(a) 


which is equivalent to 

J (Gz, 
on account of Eq. (6.2). This, together with numerical evaluations of J;, I, gives the first line of 
Eq. (6.4). 


It remains to explain the legitimacy of neglecting the term in a*£/4 in Eq. (A, 2). This term con- 
tributes to I(a) the amount 
—et 


The er converges near §=0. The [ ] has the fie: 1 for =0, and its value decreases to about 
4 for large &. The contribution is, therefore, less than 


2N 


which vanishes as Na’. The explanation of the first line of Eq. (6.4) has now been completed. The 
second line is obtained by very similar reasoning. 
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The free fall of atoms is observed in long molecular beams of potassium and cesium atoms. 

The measurement of the intensity distribution in a beam deflected by gravity represents the 

velocity distribution of the beam atoms and permits an accurate determination of this distri- 

bution. The results show that the measured values agree in general with those calculated for 

a modified Maxwellian velocity distribution in the beam. At larger deflections, i.e., for slow 

atoms, a deficiency of intensity was observed, which increased both with increasing deflection 

. and with increasing pressure in the oven where the beam originates. This deficiency is explained 
by collisions in the immediate vicinity of the oven slit. 


1. INTRODUCTION 


METHOD for the exact measure ment of the 
Bohr magneton devised recently by one of 

us! is based on the compensation of the accelera- 
tion of gravity acting on a molecular beam of 
cesium atoms by the magnetic acceleration 
produced by an inhomogeneous magnetic field. 
The effect of gravity on atoms, that is their free 
fall, is easily observable in a long molecular 
beam. In a cesium beam 2 meters long, the fall 
distance of atoms having the most probable 
velocity corresponding to an oven temperature 
of about 450°K is approximately 0.2 mm. The 
magnetic inhomogenity required to compensate 
the acceleration of gravity is in this case only 
about 25 gauss/cm. It is, therefore, possible to 
produce the magnetic field by a current in a 
straight conductor of circular cross section; thus 
the magnetic inhomogeneity does not need to 
be measured, but can be calculated from the 
geometrical dimensions and the electric current. 
Hence, it is only necessary to measure the current 
I, through the conductor which produces a mag- 
netic field of the inhomogeneity (0H/dr)o exactly 
compensating the acceleration of gravity for the 
Cs atoms with the largest magnetic moment yo 


ls 
| 
Fic. 1. Principle of the method. 
*At on leave to the National Labora- 
cago, 


in Berkeley. Ca 
Ow in e 
10. Stern, Ph 51, (1937). 


equal to one Bohr magneton. Since the com- 
pensation of the acceleration is independent of 
the velocity, this method is in principle inde- 
pendent of the velocity distribution in the beam. 
For the actual measurement of Jo, however, the 


exact knowledge of the velocity distribution in’ 


the beam is necessary. 

If the current J through the conductor is ale 
than Jo, all the atoms in the beam are deflected 
downwards. If I> Jo, those with the largest mag- 
netic moment are deflected upwards. The com- 
pensating current J») can be determined by 
placing a detector slightly above the plane of the 
beam (position D’ in Fig. 1), and by increasing 
I gradually until atoms deflected upwards strike 
the detector. In practice, the beam intensity i 
striking the detector must be measured as a 


‘function of the current J, and J is found by an 


extrapolation, which requires the knowledge of 
the actual velocity distribution in the beam. 
From published measurements by Cohen and 
Ellett? of the velocity distribution in beams of 
alkali-metal atoms undergoing magnetic de- 
flection, it appeared that in such beams the 
modified Maxwellian distribution was correct to 
within a few percent down to velocities of one- 
half of the most probable velocity. Experiments 
for the measurement of the Bohr magneton 
carried out by us from 1938 to 1942 showed, 
however, that considerable deficiencies of slow 
molecules (up to 50 percent of the number ex- 
pected on the basis of Maxwell’s law) are present 
even at oven pressures lower than those used in 


2 V. W. Cohen and A. Ellett, Phys. Rev. 52, 502 (1937). 
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the magnetic deflection experiments of Cohen 
and Ellett.* It was, therefore, necessary to deter- 
mine these deviations from the modified Max- 
wellian distribution as a preliminary work. Since 
the apparatus designed for the measurement of 
the Bohr magneton allowed such measurements 
by means of deflection by gravity (“free fall’), 
this method was used. 


2. PRINCIPLE OF THE METHOD 


To observe the free fall of molecules, the fol- 
lowing arrangement was used (Fig. 1). A beam 
of Cs atoms produced by the oven slit O and the 
collimating slit S was detected by a hot tungsten 
wire D. Both slits as well as the detector wire 
were horizontal. The Cs atoms striking the 
detector wire became ionized and the ions were 
collected on a negatively charged cylinder sur- 
rounding the wire (Langmuir-Taylor method).* 
The ion current between the detector wire and 
the collecting cylinder gave directly the number 
per second of Cs atoms impinging upon the wire. 
The dotted lines show the path of Cs atoms with 
different velocities. 

We assume OS=SD=l. Then for an atom 
with the velocity v, the fall distance is S=g/?/v* 
(not g?/2v*). For Cs atoms of the most probable 
velocity a the fall distance S, is 0.174 mmf for an 
oven temperature T=450°K and for /=100 cm. 
Since the width of the slits defining the beam and 
the thickness of the detector wire in our experi- 
ments was only 0.02 mm, this fall distance is 
easily observable. For potassium atoms, S, is 


approximately 0.04 mm, which is still well ob- 


servable under our experimental conditions. The 
intensity distribution in the beam deflected by 
gtavity measured in the vertical direction pic- 
tures the velocity distribution of the atoms in 
the beam. 

With slits of the width }, the intensity dis- 
tribution in the undeflected beam is given by a 
trapezoid (Fig. 2). Assuming Maxwellian dis- 
tribution in the oven, the fraction of atoms in 
the beam in the velocity interval between v and 


* Because of a different slit form, their “effective” 
pressure at the slit mi have been much lower than 
measured pressure. A comparison of their results 
with our results is, therefore, not possible. 
*J. B. Taylor, Zeits. f. Physik. 57, 242 (1929). 
t The value 0.177 mm given in reference 1 is incorrect. 


Fic. 2, Intensity distribution in a  idiinidlie beam deflected 
by gravity. 


The intensity at the displacement S is given by: 


OLS, 


ex 
S—2b 


S—3b Sa 
+ 

S. S—3b 
This is the intensity distribution as function of 
the displacements which would be measured with 
an infinitely thin detector wire. The finite width 
of the detector can easily be taken into account 
by integration; in our Cs experiments, this cor- 
rection amounted to less than 1} percent near 
the maximum of the intensity distribution and 
was negligible elsewhere. 

To compensate for the force of gravity by a 
magnetic force, one must produce a magnetic 
force on the Cs atoms of the same magnitude, 
but of direction opposite to the gravitational 
force. This can be accomplished by an electric 
current J underneath and parallel to the beam. 
At a distance d from the current, the magnetic 
field strength is H=2I/d, the inhomogeneity 
0H/dr =2I/d?, and the force F on an atom whose 
component of magnetic moment in the direction 
of the field is uy is F=yq-2I/d*. Compensation 
takes place if 

2Io/d* = mg, (3) 
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Fic. 3. General arrangement of the apparatus. 


where m is the mass of the atom and g the 
acceleration of gravity. Because of the nuclear 
spin (7/2) of the Cs atom and the use of ‘“‘weak”’ 
fields (not strong enough to uncouple the nuclear 
and electronic moments) the cesium beam is 
split by space quantization into 16 beams, each 
_ beam consisting of atoms with different values 
of ux. One-half of these values of uy result in a 
magnetic force in the same direction as gravity; 
for the other half, compensation is possible and 
will take place for those with the largest value 
of ug at the smallest value of 7. (One-sixteenth 
of all the atoms will have the largest value of 
#a =us—pmw, Where up is the Bohr magneton and 
un the magnetic moment of the nucleus, which 
is of the order of magnitude of 1/1000 of the 
Bohr magneton.) 

To measure Jo, the detector is placed a small 
distance above the position of the undeflected 
beam (position D’). As long as I< Jp, all atoms 
are deflected downwards and no atoms strike the 
detector. As soon as J> Jo, all the atoms with the 
largest magnetic moment (75 of the total 
number) will be deflected upwards and those in 
the proper velocity interval will strike the de- 
tector. As stated before, this method of measur- 
ing of Jo is in principle independent of the ve- 
locity distribution in the beam. However, for 
practical reasons (small intensity of deflected 
_ atoms and background of scattered atoms) the 
intensity of the deflected atoms has to be 
measured at several detector positions as a func- 
tion of J, and Jp must be found by extrapolation. 
The extrapolation requires the knowledge of the 


actual velocity distribution in the beam. This is | 


the reason for the experiments described in this 
paper. 

To find the value of uz, the quantities m, g, and 
_d in addition to J) have to be known. The mass* 


* If Eq. (3) is —e by A dro’s number N, it 
becomes Nyus2I)/d*=Mzg-2I/d*= Mg, where M is the 
atomic weight. Therefore, we measure actually the Bohr 


magneton per mole, Mz, with the accuracy discussed here. 


Fic. 4. Glass plates carrying slit system. 


m is known to an accuracy of about one part in 
10‘, and g, of course, at least to one part in 108, 
However, d is not only difficult to measure, but 
it is not simple to design an arrangement in 
which d is sufficiently defined and constant. 

If we want the result to be exact to one part 
in 1000, the distance d between the center line of 
the current and the center of the beam has to be 
exact to one part in 2000. In our apparatus, this 
distance d was 2 cm. It is not possible to make 
this distance much larger, because then J would 
be too large. Even with d=2 cm, J» is about 500 
amp. Therefore, d has to be exact to 1/100 mm. 
Even with a solid rod as conductor it is not easy 
to have the radius exact to 1/100 mm over a 
length of more than 2 meters, nor is it easy to 
hold this rod in such a fashion that the center 
line is straight and parallel to the beam within 
this accuracy. Moreover, because of the heat 
developed by the current, it is necessary to 
provide cooling. This requires the use of a tube 
instead of a rod. Even in the best precision tubes 
obtainable, the wall thickness was not uniform 
to 1/100 mm. 

To overcome this difficulty, we used the fol- 
lowing arrangement: In addition to the beam 
above the conductor, a second beam was ar- 
ranged at the same distance below. If we deter- 
mine the J» for both beams and take the average, 
this value will be the Jo for d equal to one-half of 
the distance between the two beams within the 
following accuracy: If 6 is the deviation from 
the 2-cm distance between the beam and the 
center line of the conductor, then the error in the 
average J, for both beams is of the order of mag- 
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nitude (5/d)?. The reason for this js that if the 
distance in one place is too small for one of the 
beams by the amount 4, it is too large by the 
same amount for the other beam. The distance 
between the two beams can be defined and 
measured very accurately by using a slit system 
in which both the source slits and the collimating 
slits are so designed that the edges of a single 
piece of metal form the lower jaw of the upper 


slit and the upper jaw of the lower slit. Thedis-. 
- tance between these edges was measured with a 


micrometer slide to a few thousandths of a mm. 
Even if the center line of the conductor should 
be misaligned by as much as 0.1 mm, the error 
produced hereby is only of the order of mag- 
nitude 10~*. In our experiments, the error in 
alignment was less than 0.05 mm. 


3. APPARATUS 


(a) General Arrangement (Fig. 3) 


The molecular beam was produced by Cs 
vapor streaming out from an oven slit and was 
defined by a foreslit S; and a collimating slit S». 
The introduction of a foreslit was necessary for 
two main reasons: First, in order to get a good 
vacuum in the “beam room,” it was necessary 
to separate the latter from the “‘oven room.” By 
this separation, it was possible to maintain a very 
good vacuum (10-7 mm Hg or better) in the 
beam room, where, because the beam was 2 m 
long, a long mean free path was required. In the 
oven room, the pressure was usually around 10~-° 


‘mm, but the beam length there was only 10 cm. 


Secondly, it was easier to maintain the correct 
mechanical alignment between foreslit, col- 
limating slit, and detector than between the oven 


slit and the other slits, since the oven position is 
not exactly fixed and changes during the heating. 
The alignment was the most important and dif- 
ficult part of the experiments. For the magnetic 
deflection experiments, the detector was placed 
only a few hundredths of a mm above the “unde- 
flected’ beam, since the intensity of deflected 
atoms decreases rapidly with increasing distance. 
In order to obtain the compensating current Jo, 
this intensity i was measured as a function of the 
current J producing the magnetic field. During 
these measurements, the relative position of the 
slits and the detector had to remain constant to 
about one-thousandth of a mm. Since the 
distance between the two slits and between the 
collimating slit and the detector was one meter 
in each case, the demands on the rigidity of the 
apparatus were considerable. In the first appa- 
ratus, the detector and the slits were mounted on 
a self-contained metal framework which was 
supported inside the evacuated beam room, but 
this construction did not produce the necessary 
rigidity. We finally chose the following arrange- 
ment. 

The slits and detectors were mounted on a 
slab of plate glass (‘vertical glass plate’”) 204 cm 
long, 9.5 cm high, and 1.9 cm thick, as shown in 
Fig. 4. The slits were mounted horizontally. The 
high degree of rigidity is, therefore, required only 
in the vertical direction, where the glass plate 
was very rigid with respect to bending because 
of its height of 10 cm. A slight displacement of 
the slits in the horizontal direction (up to 1/10 
mm) does not affect the accuracy. The vertical 
glass plate was resting on another glass plate 210 
cm long, 3 cm high, and 11.5 cm wide, which was 
mounted horizontally. This glass plate carried a 
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Fic. 5. (a) End view, (b) side view of apparatus. 
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Fic. 6. (a) End view, (b) front view of slit system. 


copper tube (0.947 cm O.D. and 0.47 cm I.D.) 


which served as conductor for the current pro- 
ducing the magnetic field. The horizontal glass 
plate was supported by two brass rings inside the 


vacuum envelope. 


(b) Vacuum Envelope 


The envelope consisted of a glass tube ap- 
proximately 9” in diameter with two brass end 
pieces as shown in Fig. 5. The glass tube had 
two side arms in the center for the shutter and 
for gas inlets, while the brass end pieces carried 
pump connections, traps, etc. It was divided by 
a brass plate into an oven room about 20 cm 
long, located inside one of the brass end pieces, 
and a beam room about 210 cm long. The 
oven room was evacuated by one large octoil 
diffusion pump,‘ and the. beam room by two 


similar pumps. Large liquid-air cooled traps were 


_ arranged between pumps and apparatus, The 


dividing plate had a central hole for the copper 
tube conductor and two narrow channels for the 
two beams. The central hole was equipped with 
a bushing mounted on a sylphon bellows which 
permitted a slight motion of the tube but fur. 
nished a seal of high flow resistance. The channels 
for the beams were 20 mm long, 5 mm wide, and 
1 mm high. The total flow between the two 
parts was about 1 liter/sec. as compared with 
a pumping speed of more than 50 liters/sec. for 
the beam room; that means that a pressure 
of 10-* mm Hg in the oven room produced 
only 2X10-' mm in the beam room. The 
ends were closed with plate glass disks of 26 cm 
and 1.9 cm (}”) thickness which were equipped 
with central holes of 2 cm and with bushings 
carried by sylphon bellows allowing a slight 
motion of the copper tube. All vacuum seals 
were made with Apiezon sealing compound Q. 
In order to obtain the required vacuum of the 
order of magnitude of 10-* mm, it was necessary 
to pump for several days. For the efficient con- 
densation of vapors, a copper strip 4.5 cm high, 
0.5 cm thick, and 185 cm long was suspended 
parallel to the beam from two metal Dewar 
vessels filled with liquid nitrogen. For the mag- 
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PRESSURE CONNECTION 


Fic. 8. (a) Side view, (b) end view of detector assembly and mounting. 


netic deflection measurements, this copper strip 
could not be cooled, since it did not reach tem- 
perature ,equilibrium for several hours and 
changed the position of the beam during this 
period, apparently by changing the temperature 
of the slit system unevenly. 


(c) Slit and Receiver System 


The two foreslits and two collimating slits 
were each 3 mm long and the distance between 
the upper and lower slits was 4 cm. The upper 
slits were 0.022 mm and the lower slits 0.026 mm 
wide. They consisted of a base plate A and a slit 
plate B (Fig. 6). Plate A was made of brass and 
was screwed to a brass angle which in turn was 


_ fastened to the vertical glass plate by means of 


screws inserted through holes in the glass plate. 
The plates A and B had cut-outs for the copper 
tube T and plate A had two 3-mm holes under- 
neath the slits, defining the “length” of the 
beams. The centers of the holes were 4 cm apart 
and in the same plane as the center of the copper 
tube. The upper jaw of the lower slit and the 
lower jaw of the upper slit were formed by the 
sharp edges of the center piece B. The distances 
between these edges were 3.994 cm for the fore- 
slits and 3.991 cm for the collimating slits. They 


were measured with a microscope and a Gaertner 
micrometer slide to a few thousandths of a mm. 
The upper jaw of the upper slit and the lower 
jaw of the lower slit were formed by two small 
brass pieces screwed to the base plate A. The 
width of the slits was determined by two 
methods: First, by a direct measurement with 
microscope and eyepiece micrometer, and sec- 
ondly, by measuring the distance of diffraction 
fringes with sodium light. 


(d) Detector 


The beam atoms were detected by the Lang- 
muir-Taylor method, in which each Cs atom 
striking a hot tungsten wire is re-evaporated as 
an ion and collected on a negatively charged 
plate. The plate current gives directly the 
number per second of atoms striking the wire 
surface, or the intensity of the beam. This current 
was measured with a FP-54 DC amplifier circuit 
and a galvanometer (Fig. 7). The FP-54 tube 
and the grid leak resistor of 10" were mounted 
inside the vacuum in a side tube (see Fig. 5) and 
the connecting wire between the collector plates 
and the grid cap of the FP-54 tube was self- 
supporting. Each detector consisted of a tungsten 
wire of 0.02-mm diameter and 15-mm length 
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which was mounted between springs. The col- 
lector plates consisted of nickel cylinders with 
windows for letting the beam through and were 
insulated with quartz rods. The detectors had 
to be movable in the vertical direction by about 
1.5 mm in order to measure the beam intensity 


as a function of the displacement. A movement > 


which was insensitive against vibration, free from 
friction, and reliably reproducible was obtained 
in the following way (Fig. 8). Each detector 
assembly was attached to the tip of a Bourdon 


- tube from a pressure gauge. This Bourdon tube 


was connected by a very flexible metal tube, 
wound in the form of a spiral, to a thin copper 
tube which was soldered through the vacuum 
mantle and led to a small glass bulb in which 
the pressure could be changed from vacuum to 
about two atmospheres. Calibration with a 
micrometer showed that a pressure change of 1 
cm Hg produced a displacement of the detector 
wire by almost exactly 1/100 of a mm. The 
pressure was measured with a mercury manom- 
eter and adjusted with a mercury leveling bulb 
to about 1/10 mm. The position of the detector 
wire was therefore reproducible with an accuracy 
of 1/1000 mm, the actual displacements were 
measured to better than one percent. 


(e) Alignment 


The parallelism of foreslit, collimating slit, and 
detector wire was insured by making all of them 
horizontal. Each slit plate B was provided with 
a small level bulb LZ. Both edges of the slit plates 
were machined to be accurately parallel. Then a 
surface plate was made horizontal by means of 
a sensitive level. The slit plate was put on the 
surface plate so that the edges were horizontal, 
and finally the level bulbs were so adjusted that 
the bubble was in the center position. After 
mounting the slit plates on the vertical glass 


plate, it was only necessary to bring the bubbles 
back to this position in order to guarantee 
parallelism of the slits with an accuracy of one 
part in 1000. The detector wires were made 
horizontal with the aid of a microscope with a 
cross hair, which was adjusted horizontally with 
respect to the same level. 


(f) Conductor Tube 


The copper tube conducting the current for the 
magnetic field was mounted on the horizontal 
glass plate. It rested on ten ‘‘chairs’’ spaced 29 
cm apart, which were fastened to the glass plate 
by brass clamps. The chairs were lined up by 
means of stretched tungsten wires. The tube was 
loosely tied to the chairs with thread so that it 
could move parallel to its axis to allow for thermal 
expansion. 

The position of the slits with respect to the 


.tube, that is of the two molecular beams, had 


to be invariable to a few thousands of a mm, but 
did not have to be known to better than one- 
tenth of a mm because of the two-beam arrange- 
ment. For the slit alignment, the vertical glass 
plate rested on two supports on the horizontal 
glass plate, and its position was secured with 
stops and springs. The slits were then adjusted 
so that the center lines of the beams were directly 
above and below and equally distant from the 
center line of the copper tube. This adjustment 


was made with the aid of a microscope and ad- . 


justable parallel blocks and had to be accurate 
only to one-tenth of a mm. The whole alignment 
was made outside the vacuum envelope. The two 
glass plates could be separated and reassembled 
without destroying the alignment. 


TABLE I. Ratio of measured to calculated intensities of 
the cesium beam. 


°K mm Hg Sa 3Sa Beam 
439 1.9107 1.01 0.87 Lower 
443 2.27 0.90 0.79 0.66 Lower 
449 2.95 0.93 0.75 0.67 Upper 
453 3.52 0.91 0.75 Lower 
453 3.52 0.93 0.75: 0.56 Lower 
456 3.99 0.93 0.77 0.58 Lower 
457 4.17 0.94 0.70 Lower 
476 8.87 0.79 0.51 0.33 Upper 
477 9.23 0.86 0.52 0.35 Lower 
478 9.59 0.82 0.52 0.39 Upper 
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INTENSITY 


20 30 25a 40 
Fic. 10. Gravity deflection of a cesium beam. 


(g) Ovens 


The ovens were made of monel metal as shown 
in Fig. 9. The cesium was contained in sealed-off 
glass capsules holding about 4 g. They were 
broken inside the ovens in a nitrogen atmosphere, 
and the plugs were closed immediately there- 
after.. The ovens were mounted in a cylindrical 
copper tube heated by a Nichrome coil. Another 
small heater unit was screwed on to the front of 
the oven near the slit. A thermocouple was 
attached to the back of each oven. In a separate 
test, a second thermocouple was attached to the 
front of the oven near the slit and the tem- 
perature difference between the two thermo- 


couples (about 5°C) was recorded for different 


temperatures. The ovens were mounted in a 
cradle which could be moved inside the vacuum 
in a vertical direction by means of a micrometer 
screw and a sylphon bellows. The cradle was also 
equipped with a level in order to make the oven 
slits horizontal. The size of the slits was 3 mm 
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s INTENS/TY 


WIDTH OF DETECTOR 


Fic. 11. Gravity deflection of a potassium beam. 


X0.06 mm. Each oven rested on 3 pins which 
insured the reproducibility of its position. 


(h) Assembly 


After the alignment was completed, the hori- 
zontal glass plate carrying the copper tube was 
inserted in the vacuum apparatus. Then, the 
brass plate separating the beam room from the 
oven room wag attached to the front end of this 
glass plate. The brass plate was equipped with a 
sylphon bellows and a bushing for the copper 
tube and with two channels for the beams. (See 
Fig. 5.) The cradles for the ovens were also 
attached to it. Next, the vertical glass plate 
carrying slits and detectors was slid into the 
vacuum envelope and moved into its correct 
position where it was held by two springs. Then 
the electrical connections between the detectors 
and the amplifier tube, etc. and the vacuum and 
pressure connections between the Bourdon tubes 
and the outside were made. A shutter operated 
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Fic. 


by pressure changes in a sylphon bellows was 
inserted through one of the side arms of the large 
glass tube near the collimating slit. Then the 
glass end plates provided with bellows for the 
copper tube were attached and tightened with 
Apiezon. Finally, the ovens were charged with 
Cs, inserted into the cradle through a side hole 
in the oven room and after the heater and ther- 
mocouple connections were made and the hole 
closed with a glass plate, the apparatus was 
evacuated as fast as possible. 


4. PROCEDURE 


After the necessary vacuum of 10-7? mm Hg 
or better was obtained, the oven was heated 
slowly to about 450°K and allowed to stabilize 
at a certain temperature in that region. The 
detector was moved to the position of maximum 
intensity and the oven was shifted up, and down 
until the oven slit covered the collimating slit 
completely. The beam intensity was measured 
during this operation and the resulting ‘‘oven 
displacement curve’’ served as a check for the 
slit alignment. With the oven slit in the correct 
position, the detector was moved in the vertical 
direction through the beam in steps of one or 
several hundreds of a millimeter and the gravity 
deflection curve was obtained. For deflections of 
more than 0.1 mm from the maximum, it was 
necessary to move the oven slit slightly down- 
wards with respect to the foreslit in order to 


keep all the slits, as well as the detector, on the 


parabolic trajectory of the atoms with the 
velocity corresponding to the detector position. 
The intensity of the beam was checked frequently 
by returning the detector to the maximum 
position, and a possible shifting of the position 
of the beam was controlled by checking a few 
points on the steep, upper side of the gravity 
deflection curve. The position of the “unde- 


‘12. 


flected’”’ beam was calculated from the measured 
points on the upper side of the gravity deflection 
curve, essentially from the position of the point 

at which the intensity is one-half of the maximum 
intensity. This distance S, from the upper edge 
of the undeflected beam (Fig. 2) can be calculated 
for given values of b and S,. Another method for 
the determination of the position of the unde- 
flected beam, which was also used occasionally, 
is to send a compensating current Jo through the 
copper conductor tube. The magnetic field pro- 
duced by this current will bring #5 of the atoms 
into the position of the undeflected beam, which 
can thus be measured out directly. 


5. RESULTS 


Experiments were carried out with potassium 
and cesium at different oven temperatures. The 
gravity displacement curves obtained agreed in 
general with the calculated curves. Figure 10 
shows the results for a cesium beam, Fig. 11 for 
a potassium beam. The curves are calculated, the 
dots show the measured points. For larger de- 
flections, viz., Sa, 2Sa, and 3Sa, the measured 
intensity values were definitely lower than the 
calculated values. This deficiency increased with 
increasing deflection (slower atoms) and also 
with increasing oven temperature (pressure). 


TABLE II. Values of F(y). 


F(y) —F(y* 
y=c/ce =co/c F(y) F(y)* F(y) 

0 0.5000 0.5000 

0.250 4 0.4080 0.4045 0.87 
0.333 0.3804 0.3754 1.31 
0.500 2 0.3322 0.3248 2.24 
0.667 1.5 0.2952 0.2883 

0.833 1.2 0.2737 0.2703 1.24 
0.909 1.1 0.2711 0.2697 0.51 

1 1 0.2761 0.2761 

1.1 - 0.909 0.2982 0.2967 0.51 

1.2 0.833 0.3285 0.3244 1.24 

1.5 0.667 0.4428 0.4324 2.34 

2.0 0.500 0.6645 0.6495 2.24 

3.0 1.1412 1.1261 1.31 

4.0 0.250 1.6320 1.6179 0.87 
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Table 1 shows the deficiency; expressed as the 
ratio of the measured and calculated intensities 
Imess/ cate for the positions Sa, 2Sa, and 3Sa for 
different experiments with cesium. Sa, 2Sa and 
3Sa are measured from the center of the unde- 
flected beam (not from the upper edge as in 


Fig. 2). 
6. DISCUSSION 


The experiments serve as a demonstration 
that individual atoms follow the laws of free 
fall in the same way as other pieces of matter. 
Moreover, they permit a more accurate deter- 
mination of the velocity distribution in molecular 
rays than those carried out earlier. The knowl- 
edge of this distribution is of great importance 
for many molecular beam experiments. It has 
usually been assumed that the Maxwellian dis- 
tribution law is valid as long as the mean free 
path of the molecules in the oven is several times 
as large as the width of the oven slit. These ex- 
periments show, however, that there is a con- 
siderable deficiency of slow molecule seven at 
much lower pressures. This deficiency is probably 
caused by collisions in the immediate vicinity of 
the oven slit (see appendix). 


7. APPENDIX 


(a) Calculation of the Distortion of the 
Velocity Distribution 

To estimate the order of magnitude of the distortion in 
the velocity distribution in the beam through collisions in 
the neighborhood of the oven slit (“‘cloud’’), we calculate 
this effect under the following simplifying assumptions: 
The number of collisions is calculated for rigid spheres on 
the basis of the classical kinetic theory. The cross section 
of the spheres is taken from the experimental determination 
of the mean free path. It is assumed that every collision 
throws the beam molecule out of the beam. This assump- 
tion is certainly justified in our case where a deflection as 


TABLE III. Values of F as a function of co/a. 


xo =Co/a S/Sa=xe? F 
3 1/9 0 

2 i 0.3261 
v3 0.3172 
V2.5 1/2.5 0.3148 
0.3167 
1.3 1/1.69 0.3225 
1.2 1/1.44 0.3320 
1/73 3 0.6690 
bes 
i 9 1.3319 
16 1.8827 


2 3 
Fic. 13. 


small as 10-* radian (1/100 mm to 1 m) removes the 
molecule from the beam. Furthermore, we assume that the 
density of the molecules inside the oven is uniform and 
that it and the velocity distribution of the molecules inside 
the oven are undisturbed up to the oven slit. We also 
neglect the number of molecules thrown into the beam by 
collisions outside the oven.* Finally, we consider the case 
of a circular oven opening and regard the beam, which is 
defined by the foreslit and the collimating slit, as emerging 
from the center of the oven opening perpendicular to it 
and with a cross section which is small compared to the 
area of the oven opening. 

Let / be the distance from the oven hole (see Fig. 12), 
b the radius of the oven hole, » the number of molecules 
with the velocity co passing the cross section of the beam 
at / per second, and N the number of molecules per cm? in 
the oven with the velocity c. We consider first the case 
c>co. The angle between c and co may be @. Therefore 
the relative velocity is c, = (c*+c?—2cco cos@)*. The num- 
ber of collisions at / during the time dt is: 


dn =ndt-xo*N a Cr} sindd@ (4) 


where zo* is the collision cross section and tan@y=b/l. The 
decrease of m in traversing the distance d/ = codt is therefore 


*On the basis of some rough numerical estimates we 
believe that this assumption should not change the order 
of magnitude of our results. 
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as to 20 25 
Fic. 14. 
given by: and for y=1, 
—d Inn/dt=e*(N/co) cr}sinddé (00) =(2—2 cos6o)!/6 sin%o. (7b) 
7 Limiting expressions are: 


= 
= (5) 
since 2¢,dc,=2cced cos@. Considering as function of / 
and setting we 
get 
6) 
where mo is the number of beam molecules with the velocity 
Co leaving the oven slit per second. We are allowed to 
integrate to =. because at larger distances (>b) the 
number of collisions is negligible. Substituting again for 
1, setting dl = (—b/sin*@))d@, we obtain finally: 
(6a) 
with 
S(00) = (e+e? — 2cco — (c—co)*]/6cco sin*6o. 
For the case co>c we have to replace (c—co) by (co—c). 
The numerical value of the integral 


depends only on y=c/co. For y>1, we have: 
cos6o)!— (y—1)*]/6y (7) 
for y<1 with 
cos0o)!— (y’—1)*]/6y" sin%% (7a) 
F(y)= F(y')/7’ 


1(slow beam molecules) F(y) = $y — 1/8 = }y —0.3927 
y\(fast beam molecules) = $—0.3927y, 


Table II and Fig. 13 give some values of F(y).* The 
weakening of the beam molecules with the velocity c» by 
collisions with oven molecules of the velocity c is therefore 
given by: 

F=F(y)=F(c/co), (8) 


where N is the number of oven molecules per cm* with the 
velocity c. If we take c as the average velocity and N as 
the total number No of molecules per cm’ we get already 
a fair approximation, especially for the case of slow beam 
molecules. It is, however, not difficult to take into account 
the velocity distribution of the oven molecules. We simply 
replace N according to Maxwell by 


dN= 


* F(y) was calculated by numerical nennaiion and by 
the formula 
— Cy 188, 


6F = F(x, r/2)—F(«, 
bE=E(x, /2)—E(«, «/4); 


where 


and 
e=4y/(y+1). 
Sar integrals of the first and the second 
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VELOCITY DISTRIBUTION 


Since the e-factors of the molecules with different velocities 

multiply, we get: 

=exp[—r0%NoF] (9) 

where 

4 

To evaluate F, the following approximation formula F* 


was used for F: 


For y<1;_ 
y=c/co 
This formula gives the right limiting expressions for small 
and large values of y. The factor 0.16884 of the third 
term is so chosen that the value of F for y=1 is correct. 
In the intervals 0<7y<1 and 1<7< ©, the largest devia- 
tion from the correct value is less than 24% (compare 


Table II). 
The resulting values of F as function of co/a are given 


in Table III and Fig. 14. 


(b) Comparison with the Experiments 
Our theoretical formula gives the “weakening factor” ® 

b=n/no 
The cross section for collisions of Cs atoms with Cs atoms 
was measured by Foner,:* 

=2.35 X10- cm*. 

The number N of Cs atoms per cm in the oven is: 

N=p/kT=2.1X10"Xp’, 
where p’ is the pressure of the Cs vapor in the oven in 
units of 10-? mm Hg. For the temperature we took an 
average value T= 460°K. 

The width of the oven slit in our experiments was 
5X10-* cm. In our theoretical equation, b is the radius of 
a circular oven opening. If we set b=5 X10? cm we should 
get the right order of magnitude. 

With these numerical values we have: 

(n/no) exp[—0.25p’F J. 


To compare our experimental results with this formula 
we have to consider that our Imeas/Jesic (see Table I) 


*S. N. Foner, Thesis, Carnegie Institute of Technology, 
Pittsburgh, 1945 (see also the following paper). 
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TaB.e IV. Calculated values of Jmeas/I calc. 


Sa 2Sa 3Sa 
2 0.98 0.90 0.84 
4 0.95 0.82 0.70 
8 0.90 0.67 0.50 


refers to an Jeaic which sets arbitrarily @=1 in the neigh- 
borhood of the maximum of the deflection curve. In 
reality we know only that ® is approximately constant in 
this region. This follows directly from our experimental 
result that in the neighborhood of the maximum the shape 
of the measured intensity curve agrees with the shape 
calculated from the modified Maxwellian distribution. It 
also agrees with our theoretical result for F (see Fig. 14). 
We assume therefore that our Jcaie corresponds to the 
minimum value of F=0.32. If we call the weakening 
factor in this region #9, we have: 
X0.32] =exp[—0.080p’]. 


For the deflection Sa, the velocity of the beam molecules 
is a, and from our curve for F we take PF =0.37, and 
finally obtain 


= exp[—0.25 Xp’ X0.37] =exp[—0.0925p’]. 
Hence, for the deflection Sa, our I meas/Jeate is 
T meas/Teste= exp[ — 0.01259"). 
Correspondingly we have for the deflections 2Sa and 3Sa: 


2Sa: ¢o=a/v2; F(1/v2)=0.52; 
T meas/Tcaic = exp[ —0.0509'] 


3Sa: ¢o=a/V3; F(1/v3)=0.67; 


The following Table IV gives some calculated values of 
I meas/ Teate: 

A comparison with the measured values in Table I 
shows agreement in the dependence on pressure and 
velocity as well as in the order of magnitude of the numer- 
ical values. Generally, the measured deviations are some- 
what larger than the calculated ones. That may be due to 
additional collisions, either with foreign gas molecules in 
the oven room, or directly at the oven slit with foreign 
gas molecules originating in the oven. The fluctuations 
in the measurements corroborate this explanation. The 
discrepancy could also be due, at least in part, to the 
choice of 5. A slightly larger value of (74 107*) would 
give agreement within the limits of error of the measure- 
ments. Considering the number of serious simplifications 
made in the calculations we cannot expect a better agree- 
ment than in the order of magnitude. . 
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The Mean Free Paths of Cesium Atoms in Helium, Nitrogen, and Cesium Vapor 


I. EsterRMANN, S. N. Foner,* AND O. STERN** 
Research Laboratory of Molecular Physics, Carnegie Institute of Technology, Pittsburgh, Pennsylvania 
(Received November 29, 1946) 


The mean free paths of cesium atoms in helium, nitrogen, and cesium vapor were measured 
with a molecular beam apparatus permitting the measurement of scattering deflections of 
about 5 seconds of arc. Utilizing the free fall of the beam atoms, the variation of the mean 
free path with the velocity of the beam atoms could be determined directly. This variation 
was found to be in agreement with the equations calculated on the basis of classical theory 
treating the atoms as elastic spheres. The measured values of the sums of the effective atomic 
radii are 12.0 10-* cm for Cs-He, 17.2 10-* for Cs-N2, and 27.3X10-* cm for Cs-Cs colli- 
sions. These values correspond to mean free paths of 2.1 meters in He, 2.2 meters in Nz, and 
1.36 meters in Cs-vapor of 10-*-mm Hg pressure for cesium atoms with a velocity corresponding 
to the maximum intensity of the gravity deflection curve. ‘ 


S pointed out in a previous paper,’ the 

deficiency of slow atoms in a molecular beam 
of Cs-atoms can be explained as arising from the 
scattering of the beam atoms in the immediate 
vicinity of the oven slit. For a quantitative 
estimate of the effect of this scattering, it is 
necessary to know the effective cross section for 
the Cs-Cs collision. The present investigation 
was undertaken to measure the mean free path 
of cesium atoms in cesium vapor, helium, and 
nitrogen, and to determine the velocity de- 
pendence of the mean free path. ' 


atoms with the DeBroglie wave-length \=h/ Mp 
incident upon an atom which is at rest in the 
coordinate system. The scattering cross section 
Q is defined as 


f I(@) sinédé, (1) 
0 


where J(@) is the scattering intensity per unit 
solid angle of the beam atoms scattered through 
an angle @ from the direction of the incident 
beam. For a hard sphere model, Massey and 
Mohr? have given as a rough approximation the 


* Now at footed 
** Now in keley 


1. THEORY 


In the quantum-mechanical treatment of scat- curve A. Here ro is the sum of the “classical” 
tering, Massey and Mohr* consider a beam of atomic radii and k=2x/A, where \=h/ Mog is 


Fic. 1. Scattering intensity 
as function of angle. 


ty Laboratory, Johns Hopkins University, Silver Spring, Maryland. 
, California. 

11, Estermann, O. C. Simpson, and O. Stern, Phys. Rev. 71, 238 83k referred to as I. 
*H. S. W. Massey and C. B O. Mohr, Proc. Roy. Soc. 141, 454 (1933 
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angular dependence of J(@) shown in Fig. 1 as 


1947 
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Fic. 2. Schematic diagram of apparatus. 


the DeBroglie wave-length associated with the 
system; M =m m2/(m,+mz2) being the reduced 
mass and vg the relative velocity of the atoms. 
The large increase in scattering intensity for 
small angles makes the cross section Q about 
twice the classical value ro*. Curve B in Fig. 1 
represents J(6) siné. The experimentally meas- 
ured cross section depends on the resolving 
power of the apparatus. If a scattering corre- 
sponding to an angle 6’ can be detected, the 
error made in measuring the total scattering 
cross section is found by comparing the unshaded 
area under curve B in Fig. 1 with the total area. 
The resolving power necessary to get meaningful 
results is reached if an increase in the resolving 
power will not materially increase the measured 
cross section. In our experiments, this condition 
was fulfilled. 

If an attractive potential of the form V= 
—C/r§ exists, the total cross section is given 
by? 

Q=[(25—3) (2) 


where k is again 27/\ = 24 Movpg/h and f is a known 
function of s. Assuming dipole-dipole Van der 
Waals forces with a potential V= —C/r®, Q takes 
the form 

Q= Bor", (3) 


where B is a constant, so that the cross section 
varies as the inverse 2? power of the relative 
velocity. The expected dependence of Q on the 
relative velocity is therefore a slowly varying 
function, with the result that the variation of Q 
with the beam atom velocity is still smaller, so 
that for the range of velocities encountered in 
the experiments we can consider Q a constant, 
especially in view of the accuracy of the measure- 
ments and the degree of monochromization of 
the beam. 


*H. S. W. Massey and C. B. O. Mohr, Proc. . Soc. 
144, 188 (1934). 
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Fic. 3. Scattering chamber for nitrogen and helium. 


2. EXPERIMENTAL ARRANGEMENT 


The mean free path was determined by meas- 
uring the weakening of a molecular beam of 
cesium atoms passing through a scattering 
chamber filled with helium, nitrogen, or cesium 
vapor. The apparatus used for these experiments 
was the one described in paper I, with a scat- 
tering chamber C mounted between the foreslit 
and the collimating slit of the upper beam, as 


‘shown in Fig. 2. For helium and nitrogen as 


scattering gases, the scattering chamber was 
made of two brass pieces (Fig. 3). The lower 
piece had a milled groove 0.4 mm deep and 5mm — 
wide, which together with the upper piece 
formed two channels 0.4 mm high, 5 mm wide, 
and 1 cm long through which the beam passed. 
The large flow resistance of these channels and 
the high pumping speed in the rest of. the ap- 
paratus permitted the maintenance of a pressure 
ratio of about 1000 to 1 between the chamiber and 
the apparatus. The gas was fed into the chamber 
through a flexible pipe line from a 5-liter flask 
outside the apparatus. The pressure in the 
chamber was controlled by a sensitive needle 
valve and measured with an ionization gauge 
connected to the chamber through another tube. 
A by-pass permitted the outgassing of the 
ionization gauge. 

For the cesium-cesium experiments, the scat- 
tering chamber had to be made of a metal not 
affected by cesium. Monel metal was found to 
be satisfactory. As shown in Fig. 4, this scattering 
chamber was machined from a solid metal block. 
The scattering took place in a tube T 3.5 cm 
long and 0.476 cm in diameter through which 
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HEATER 


Fic. 4. Scattering chamber for cesium vapor. 


the beam passed. The cesium vapor was supplied 
from a cesium capsule in the second part of the 
scattering chamber, which was connected to the 
first by two }-inch holes. Since the cesium vapor 
escaping from the scattering chamber was 
quickly condensed at the cooled copper strip (see 
Fig. 5 of paper I), it was possible to use slits S; 
and S; instead of channels at the ends of the scat- 
tering tube without materially increasing the 
pressure in the apparatus. A Nichrome ribbon 
heater was fastened to the top of the chamber 
and a thermocouple to the plug. For better 
thermal insulation, the chamber was mounted 
on a glass plate and supported by brass angles 
attached to the vertical glass plate (see Fig. 4 
of paper I), which also carried the slit and 
detector system. 

It was found that the Cs pressure in the scat- 
tering chamber could not be calculated reliably 
from the temperature measured by the thermo- 
couple. In order to measure the pressure of the 
scattering Cs vapor, a third slit Ss was attached 
to the side of the scattering tube and am auxiliary 
hot tungsten wire detector was mounted 2.3 cm 
away from the side slit. A collimating slit and a 
magnetically operated shutter were mounted 


between this detector and the side slit. From the 
intensity of the cesium beam striking the side 
detector and the geometry of the arrangement 
the pressure of Cs atoms in the scattering tube 
could be calculated. It was also calculated from 
the intensity of the Cs beam which originated 
from the scattering chamber. This intensity was 
measured with the beam detector after displacing 
the oven sufficiently as to cut out all the atoms 
originating in the oven. Because of the efflux of 
cesium vapor through the slits at the ends of the 
scattering tube, the pressure inside the chamber 
was not quite uniform. An appropriate correction 
for this effect was applied in the calculations. 


3. RESOLVING POWER 


The resolving power of the apparatus follows 
from the geometry of the molecular beam. The 
source slit (fore slit) and collimating slit were 
both 0.02 mm wide and 100 cm apart. The 
detecting wire had a diameter of 0.02 mm and 
was located 100 cm beyond the collimating slit. 
The distance between the middle of the scat- 
tering chamber and the collimating slit was 85 
cm. Hence, if the direction of the velocity vector 
of a beam atom passing through the scattering 
chamber is changed by an angle of 2.410 
radian or 5 seconds of arc (0.02 mm in 85 cm), the 
atom cannot pass through the collimating slit 
and is, therefore, measured as “‘scattered.’’ The 
actual resolution will be slightly lower since 
atoms can be scattered into the beam from the 
0.1-mm broad beam wedge inside the scattering 
chamber, which is formed by the oven slit and 
the fore slit. 

In the theory of scattering, the angle of scat- 
tering is the angle through which the relative 
velocity vector vg is turned. What is measured 
experimentally is the deviation of the beam 
atom. In the case of small angular deviations, 
this is roughly half the angle through which the 
relative velocity vector is turned.‘ 

By using the hard sphere model and the 
Massey and Mohr approximation (Fig. 1) for 
Cs-Cs scattering, the angle 6)=2/kro at which 
quantum scattering begins to make itself 
noticed, turns out to be 0.3°, if the effective 
radius for the Cs atom is taken from our experi- 


‘W. H. Mais, Phys. Rev. 45, 773 (1934). 
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ments as equal to approximately 10X10-* cm. 
This means that beam atom deflections of less 
than 0.15° will show “quantum scattering”’ (see 
Fig: 1). Since the apparatus has a resolving 

100 times larger than this, the interpre- 
tation of the measurements seems justified, even 
allowing for the fact that there is no a priori 
reason for assuming that large atoms like Cs will 
behave in collisions like hard spheres, and that, 
therefore, the actual form of the angular function 
](6) is not known. 


4. PROCEDURE 


(a) Helium and Nitrogen Experiments 


With the scattering chamber evacuated, the 
beam was stabilized in the manner described in 
paper I, and the gravity deflection curve of the 
unscattered beam was measured out. Then, a 
certain gas pressure was introduced into the scat- 
tering chamber. The measurement of the gravity 


"curve was repeated for different scattering 


pressures between 8 and 33 X10-* mm Hg. 
(b) Cesium Experiments 


For the cesium experiments, the scattering 
chamber was charged with a Cs capsule which 
was broken in a Nz atmosphere before assembling 
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the apparatus. After a vacuum of better than 
10-* mm was obtained, both the oven and the 
scattering chamber were heated slowly to a tem- 
perature somewhat higher than the operating 
temperatures and outgassed for several hours. 
Then the scattering chamber was allowed to 
cool off to room temperature and the unscattered 
beam was measured out. The scattering chamber 
was then heated again until the beam intensity 
was cut to $ or } of its original value; and after 
stabilization, the gravity curve was measured 
again. Since part of the intensity measured by 
the beam detector came from cesium atoms 
originating in the scattering chamber, the oven 
was later displaced sufficiently as to cut out all 
the atoms originating from the oven, and the 
residual intensity was measured out. These 
measurements served as a second determination 
of the pressure in the scattering chamber. In 
some experiments, a check of the original beam 
intensity was made ~by letting the scattering 
chamber cool off again to room temperature. 


5. RESULTS 
(a) Helium and Nitrogen Scattering 


Gravity deflection curves obtained in the 
helium and nitrogen experiments are shown in 
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‘Fic. 5. Intensity distribution in cesium beam scattered by helium 
(o measured points, x calculated points). 
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BEAM INTENSITY CARBITRARY UNITS) 
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Fic. 6. Intensity distribution in cesium beam scattered rm nitrogen 
(o measured points, x phe points). 


Figs. 5 and 6, respectively. In both figures the 
higher intensity curve represents the “‘unscat- 
tered’”’ beam. Similar curves were obtained with 
helium pressures of 8.5 X 10-* mm and 
mm Hg, and a nitrogen pressure of 18.010-° 
mm Hg. 

The points marked ‘‘calculated” were found 


in the following way: The weakening of a. 


homogeneous beam (i.e., a beam in which all 
atoms have the same velocity) from an intensity 
I, to an intensity J is given by 
(4) 
where p is the scattering pressure, / the length 
of the scattering chamber, \ the mean free path, 
and \»=A/p a constant independent of pressure. 
If we consider the actual beam as a superposition 
of a number of homogeneous beams of different 
velocities c, we can write | 
exp(—lp/o’), (4a) 
where the superscript c denotes that the symbol 
pertains to the velocity c. If the weakening of 
the beam by a certain scattering pressure p is 
known for two velocities c and c’, then 
In(Jo*/I*) 


©) 


Because of the finite width of the slits, all the 
atoms striking the detector at a given deflection 
position S do not have the same velocity. Cal- 
culations given elsewhere® show, however, that 
the average of the free paths of all these atoms 
does not deviate appreciably from the free path 
of the atoms reaching the position S and origi- 
nating from the center of the undeflected beam. 
We may, therefore, define an effective velocity ¢ 
corresponding to every detector position S, and, 
since the deflections are inversely proportional 
to c; c/a?=S,/S, where a refers to the most 
probable velocity of the beam atoms, and S, to 
the corresponding deflection. Comparing all the 
free paths with the free path of the atoms in the 
maximum of the gravity curve, we get 
In(Io°/I*) 
doe In(Jom/I™) 
where the superscript m, refers to the atoms in 
the maximum of the gravity curve. If the rela- 
tionship between Ao* and A” is known, the scat- 
tered curve can be calculated from the unscat- 
tered curve if one —- of the scattered curve is 
known. 


N. Foner, Thess, Carnegie Institute of Technol 
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TaBLe I. Mean free paths for cesium atoms in helium 


and nitrogen. 
gas 
Heli Nitrogen 
3.0 5.70 2.38 0.268 1.460 0.346 1.212 
64 2.67 1.63 0415 1.000 0.620 1.000 
28 134 1.16 0.148 0.711 
144.1 100 1.00 0.069 0.613 0.125 
22.7 0.67 0.82 0.019 0.503 0.040 0.586 
28.2 0.50 0.71 0.007 0.435 0.017 0.519 


For helium scattering, the cesium atoms 
traversing the scattering chamber are essentially 
stationary targets for the much faster He atoms, 
hence the probability of a cesium atom being 
scattered out of the beam is simply proportional 
to the time spent in the chamber t=//c. There- 
fore, 4o"/Ao° =c"/c, and the intensity distribution 
in the scattered beam is given by 


In(Iot/I*) = (c*/c) In(Io"/I™) 
= (Sa/S")'/(c/a) In(Io"/I"). (6) 


For the collisions between cesium atoms and 
nitrogen molecules, the relationship between free 
path and velocity is more complicated. Calcula- 
tions for this case are given elsewhere.*-? 

The absolute value of the mean free path of 
the cesium atoms is given for those atoms which 
make up the maximum of the gravity curve, 
since there the measurements are most accurate. 
For these atoms, and a scattering pressure of 
10-* mm Hg, the mean free path is 


Ao” = 209 cm/10-* mm Hg for cesium in helium, 
and 
Ao” = 216 cm/10~* mm Hg for cesium in nitrogen. 


The calculated values for the intensity J* in the 
scattered beam as function of the deflection S 
are indicated on Figs. 5 and 6. They are also 
given in Table I together with the relative values 
for the mean free path for cesium atoms of 
different velocities in helium and nitrogen. 
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Fic. 7. Intensity distribution in cesium beam scattered by cesium vapor 
(0 measured points, x calculated points). 


- Meyer, Kinetic Theory of Gases (London, 1899 


E 
. H. Jeans, Dynamic Theory of Gases (Cambrides 1916) second edition. 
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(b) Cesium-Cesium Scattering at S=25S,). This effect is discussed in detail ” 
Gravity deflection curves obtained in Cs-Cs i Paper I. Otherwise, the agreement between 
scattering experiments are given in Fig. 7. Curve ™easured and calculated values is very good. 
I chawse the unwealsened beam. curves II and III For the determination of the absolute value of 
the beam passing through the cntiusiog chamber the mean free path of cesium atoms in cesium 
vapor, a series of measurements of In(J»”/J™) in 


filled with cesium vapor of a pressure of 4.9 10~* 
and 12.4xX10* mm Hg, respectively. The 
“‘undeflected’’ beam curve, as calculated from the 
slit geometry and the maximum of the “un- 
weakened” beam curve and indicated by the 
trapezoid, has been reduced by a factor of three 
to allow a better representation of the gravity 
curves. The calculated points of the weakened 
beam curves were obtained by the same method 
as in the nitrogen scattering experiments. The in- 
tensity of the unweakened curve as calculated 
on the basis of a Maxwellian velocity distribution 
in the undeflected beam indicates a deficiency of 
slow molecules (9 percent at S=S;, 30 percent 


the maximum of the gravity curve was carried 
out, with the pressure in the scattering chamber 
varying between 35 and 200X10-* mm Hg. 
These results are given in Fig. 8, while Fig. 9 
shows the resulting values of \»” as a function of 
the scattering pressure. The mean free path )y" 
at a pressure of 10-* mm Hg is 136 cm. Values for 
other velocities are given in reference (5) and can 
be calculated from Eq. (7). 


(c) Collision Cross Sections for Cesium Atoms 


From the mean free path at a given pressure, 
the classical collision cross section ro4¢* can be 


| 
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calculated’ by means of Eq. (7), 
(7) 
where 


and vg is the number of scattering atoms per 
cm’, ca the velocity of the beam atoms, T¢ the 
absolute temperature of the scattering gas, k 
Boltzmann’s constant, m4 and mg the mass of 
beam and scattering gas atoms, respectively, and 


v(x) =xe~** + (2x?+-1) f e~vdy. 
Since 
and 


where the subscripts G and A refer to the scat- 
tering gas and the beam atoms, respectively, we 
have x=(ca/as)(Tame/Tqma)}. For ca, we 
choose the velocity of the atoms in the maximum 
of the intensity curve. From this curve, we take 
¢4=1.64a4. Using the data for corre- 
sponding to this velocity, we obtain the values 
for the cross sections given in Table II, where 
cag is the sum of the effective radii, ro4¢* the 
quantum-mechanical collision cross section and 
ro the classical sum of the atomic radii calculated 
under the assumption? that the classical cross 
section is one-half of the quantum-mechanical 
cross section. 

The values of the mean free path are estimated 
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TABLE II. Effective collision radii and cross sections. 


X10") =—s re 4g* X10"* re X10* 
cm cm cm 


Encounter cm 
Cs-He 142 12.0 446 8.5 
Cs-N: 298 17.2 936 12.1 
Cs-Cs 743 27.3 2350 19.3 


to be accurate within 10 percent, hence the 
accuracy of the effective radii should be 5 per- 
cent. It is interesting to compare the classical 
radii for the Cs-Cs and Cs-He collisions. The 
radius of the Cs atom in the latter is 9.6 10-* 
cm, which is larger than the sum of the radii of 
the Cs-He collision, which is 8.5 10-* cm. This 
means, of course, that attractive forces are 
present in the Cs-Cs encounters. 

The effective collision radius for Cs-He en- 
counters was measured by Rosin and Rabi,* who 
found o4g¢=7.18X10-* cm, which is about 40 
percent lower than our value. Their apparatus 
was much smaller (beam length about 10 cm), 
and the resolving power of the order of 1 minute 
of arc as compared with 5 seconds of arc in our 
apparatus. A possible explanation for this dis- 
crepancy is that the scattering intensity for small 
angles may be much higher than is expected on 
the basis of the simplified quantum theoretical 
treatment of the atoms as hard spheres. 


*S. Rosin and I. I. Rabi, Phys. Rev. 48, 373 (1935). 
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Electron Waves in the Magnetic Dipole Field of a Neutron 


CuarLes L. CRITCHFIELD 
University of Minnesota, Minneapolis, Minnesota 
(Received October 21, 1946) 


The relativistic wave functions for an electron moving 
in the field of a central magnetic dipole that has spin $4 are 
studied. Spin inertia and retardation of the central field 
are neglected. In general, there are then four, first-order, 


coupled differential equations for the radial dependence — 


of the electron waves; and these are equivalent to two 
coupled Schroedinger equations having a positive definite 
potential energy operator. In the particular case that the 
total angular momentum of the system is zero, however, 
the equations are equivalent to a single Schroedinger 
equation in which the potential energy is positive at all 
radii only for the electron-neutron interaction. The 
positron-neutron potential becomes strongly negative in a 


small region near a separation of the order of e*/Mc*, where. 


M is the mass of the neutron. The attraction is just suf- 
ficient to form a bound Pj-state of the positron with zero 


binding energy. The reality of such neutron-positron state 
is open to question because of the simplifications that have 
been made and also because the magnetic dipole moment 
of the neutron is most probably of extended origin. The 
single equation obtained when the total angular momentum 
vanishes is the Mathieu differential equation which may 
be solved in terms of an integral function by a method 
developed by Lindemann. Asymptotic forms of these soly- 
tions, as well as of the fourth-order equation for total 
angular momentum unity are presented. The results 
indicate agreement with those obtained in various appli- 
cations of perturbation methods. In particular, the am. 
biguity of the wave functions discussed previously by 
Bloch is apparent in the solutions for unit angular mo- 
mentum. 


1. INTRODUCTION. THE WAVE EQUATIONS 


THEORY of elementary particles must 
contain, to some degree, an interpretation 

of electromagnetic fields near singularities. An 
interesting portion of the singular behavior lies 
in the magnetic field arising from the dipole 
_ carried by the particles. For électrons and posi- 
trons, and presumably for charged mesons, the 
magnetic field is concomitant with the singular 
electric field. For protons and neutrons, most of 
the magnetic field is thought of as arising from 
an extended source characteristic of the inter- 
action between nucleons and mesons. In either 
case, further knowledge of the properties of the 
magnetic fields is desirable. It is the purpose of 
this paper to present a theoretical investigation 
of a simple manifestation of magnetic dipole 
' forces, viz., the relativistic motion of an electron 
in the field of a point dipole. The line of approach 
taken here is chosen to be as parallel as possible 
to the familiar treatment of two-body problems 
where electric forces are dominant, i.e., the 
hydrogen atom and scattering by the Coulomb 
field. It is not the purpose of this work to replace, 
or refine, previous considerations made by many 
authors on hyperfine structure of spectral terms, 
magnetic scattering of electrons by nuclei and 
the scattering of neutrons by ferromagnetic sub- 


stances for which perturbation methods are 
adequate. 

If the magnetic field in which the electron 
moves were due to a magnetic dipole, w, fixed at 
the origin of coordinates, the potential acting on 
the electron would be the vector potential, A, 
which does not vary with time, 


A=(yXr)/P. (1) 


The wave equation for the electron would then 
be found by putting A in its proper place in the 
relativistic Hamiltonian function, as given by 
Dirac’s theory. Then, if the self-energy of the 
dipole field is neglected, the situation is quite 
analogous to the result of the usual procedure 
when the central force is electrostatic (over- 
looking the fact that angular momentum will 
not be conserved under the postulates made so 
far). 

If the magnetic dipole is associated with an 
elementary particle that is heavy compared with 
the electron, it is appropriate to fix its position 
at the origin of coordinates but not to fix its 
direction in space. The direction of the dipole is 
rather to be correlated with the spin axis of the 
central particle. In the present work we consider 
only the case for which the spin of the central 
particle is $4. Using the Pauli spin vector oper- 
ator @, the y in Eq. (1) must be considered as an 


258 


8 


| 
and 
the 
ore 
witl 
vari 
pote 
inte 
indi 
low 
pres 
erti 
alon 
criti 
inte 
bilit 
prot 
neut 
by t 
so tl 
state 
of tl 
and 
Han 
elect 
of gr 
(Clar 

| 


947 


ELECTRONS IN A MAGNETIC DIPOLE FIELD 259 


operator on the spin states of the heavy particle 

(2) 
where v’ is a dimensional constant. Since we also 
assume that the central particle does not carry 
an electric charge, the closest physical system to 
the one postulated here is the electron-neutron 
system. Comparison of the theory, which is for 
point dipoles, with neutron-electron systems 
suffers the obvious defect that the neutron’s 
magnetic moment is probably not of point origin 
but rather arises in a region of radius comparable 
to the range of nuclear forces. 

Correlating the central magnetic moment with 
the spin of the central particle means that the 
wave function depends not only on the position 
and spin coordinates of the electron but also on 
the spin of the neutron. It also means that the 
orientation of the magnetic field may change 
with time and this has two effects. First, the 
magnetic field must have some rotational inertia 
that should be added to the kinetic energy of the 


" system and, secondly, the vector potential now 


varies with time after the fashion of retarded 
potentials. Generally speaking, the magnetic 
interaction is so weak that angular velocities 
induced in the central field should be relatively 
low and the effects of motion negligible. In the 
present paper, consequently, the rotational in- 
ertia and retardation effects will be disregarded, 
along with the self-energy of the field. The main 
criticism of this procedure is that for the most 
interesting consequence of the work, the possi- 
bility of bound states, the effects neglected are 
probably not small. On the other hand, bound 
neutron-positron systems are probably ruled out 
by the finite extension of the neutron’s moment 
so that the whole question of the reality of bound 
states is better postponed until a study is made 
of the magnetic interaction of more elementary, 
and therefore charged, particles. 

Under the conditions just set forth, the 
Hamiltonian function that is relativistic for the 
electron! but neglects the motion of the center 
of gravity of the neutron may be given the form: 


H=c(a, Pp) +e(a, A) +Bme’, (3) 
A=—exr. 


'P. A. M. Dirac, Principles of ntum Mechanics 
(Clarendon Press, Oxford, 1935), second edition, Chap. XII. 


Neither the angular momentum of the neutron, 
3eh, nor the total angular momentum of the 
electron M=m-+4e,h commutes with H so that 
neither, alone, is a constant of the motion. 
However, it can be shown that H does commute 
with the total angular momentum 


This means that the operator H couples the 
various states that can be formed with a given 
total angular momentum. In the terminology of 
hyperfine spectroscopy, the total angular mo- 
mentum of nucleus plus electron is designated 
Fh. With a nuclear spin of $4 and a single 
electron, the angular momentum in the electron 
state, Mh, may have either of two values given 
by M=F+} and M=F-—4, if F>0; or just 
M=} when F=0. The complete wave function 
for a particular value of F will be a sum of 
product wave functions, each product containing 
one polarization of the electron wave and one 
polarization of the neutron wave. Each product 
is itself composed of the four components of the 
electron wave. Let the neutronic spin function 
be x., where s= +} designates the z-component 
of the spin state. Each product may thus be 
symbolized by 
The only term in H that operates simultane- 


ously on the y’s and x, is that containing A and 
this may be written 


sin@ cosd 
oy siné sing 
Qe Gz cosé 


More explicitly, the operation of @Q on 
v2, Ws, Val Xe may be written : 


v2, vs, Wal xe 
=isx, { —e~ es, yo, sind 
(2s+1)y2, (2s—1)y1} cosé 


1 
X (2s cos#+ sing) sin@. (4) 
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With this formula, it is a matter for direct 
computation to determine the wave equations 
for the electron. For that purpose, we adopt the 
method presented by Darwin? who expresses the 
angular dependence of the electron wave func- 
tions by writing the whole harmonic as 


Py=(—)*(k—u)! 
d k+u 
d cosé 


2*k! 


The notation differs from that in the reference 
only by the factor (—)“ which arises from a 
different choice of representation of the rotation 
group. 

Let a and 6 be functions of r alone, then one 
electron wave for which M=k+4=M, has the 


general form 
{(2k+1)P,*a, 0, 


Another form may be derived from this by 
interchanging function 1 with 3, and 2 with 4, 
and changing the sign of a relative to 6. The 
latter is equivalent to transformation by the 
matrix —ia,aya, which commutes with every 
member of H except Smc* with which it anti- 
commutes. Wave equations that one set of elec- 
tron waves obeys can then be obtained from those 
that the other obeys by substituting — mc? for mc 
in H and identifying b with a and a with —d. 
A state of total angular momentum Fh and 
z-component F,h, with F= F,=k-+1 is then the 
simple product 
b) = {(2k+1)P*a, 0, 
(5) 


The state with which 7; is coupled by the 
operator Q is of the form 


Sulf-g) =40 0}x-1 
i(2k+2) 
(6) 


where f and g are functions of r alone. It can 

be verified, that under infinitesimal rotations, 

both 7; and S, belong to the (2k+3)—dimen- 

sional representation of the group of rotations 

and that they belong to the same representation 
2C. G. Darwin, Proc. Roy. Soc. A118, 654 (1928). 


of the reflection group.* Interchanging y, with 
Vs and v2 with VW leads to two other total wave 
functions that belong to the same representation 
of the rotation group but to the alternate repre- 
sentation of the reflection group. Since Q jg 
invariant under inversions these two types of 
total wave functions are not coupled by Q. 

Operation by Q on 7; and S, can then be 
calculated by using Eq. (4) and, with F=k+1, 
leads to 


— (2F+1)QT;(a, b) 
= FT,(0, a) + Fi(F+ 1) 4S,(b, a), 


— (2F+1)QS.(f, g) 
= f)-+(F+1) Sila, f). 
Note that the order of the arguments (a, d), etc., 
is interchanged by Q. We introduce the char- 
acteristic length, w, 
w= v(e*,’ 


(v is just the number of nuclear magnetons in — 


the spin direction) and combine Eq. (7) with 
the general considerations of reference 2 as they 
relate to Eq. (3) to obtain the differential equa- 
tions in r that a, b, f and g must obey: 


da F-1 F w 

—a 

dr +r  2F+ir 
(F(F+1))! w 
—f + 

2F+1 he 

F w 

2F+1 

(FF +1)! w 


+——a=0, 


=0, 


dr r 


(8) 


+— 
F(F+1))# mec—E 
(F+1)) 
2F+1 he 
df F+iw 
dr or 
(F(F+1))! w me+E 
—at g=0, 
2F+1 he 


werk Wigner, Gruppentheorie (Vieweg, Braunschweig, 


A 


(8) 
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where E is the eigenvalue of H. Another set of 


equations for the same eigenvalue is obtained by 
changing the sign of m and exchanging the 


in Eq. (10); let 
= mic 


components of the electron wave function as then the second-order equations are: 


described above. 


Equation (8) can also be derived by a method parallel 
to that used by Dirac, reference 1, in obtaining the equa- 
tions in spherical coordinates for the hydrogen atom. This 
may be done by defining an operator 

Kh=(¢, M)+h, 
where @ still refers to the neutron spin, and making use of 
Dirac’s operator 
€=pi(@,, 


It is readily verified that 
Q=(Ke—eK)/2i, 
and from these relations, and the fact that Q is a two-rowed 
matrix in S;, Ti-space, Eq. (8) may be derived. This 
method was not used, here, as it is less clear how it may 
be extended to magnetic interaction with particles of 
The differential Eq. (8) may be reduced to a 
more compact form by the substitutions: 


y= 
(F+1)ig], 
(9) 
¢3=r(Fif—(F+1)ia], 
Fia+ (F+1)4f], 


from which one obtains the first-order, coupled 
equations for the ¢’s 


dg, fw 1) 
dr 
do: ¥ 0, 
dr 
(10) 
d¢s y 
dr fr hc 
do /1 


dr 


The equations may also be expressed as two sets 
of coupled, second-order equations which are 
readily derived by differentiation and elimination 


“) 


de 


(12) 


-1(= 


Each of these sets of two equations can be 
regarded as a pair of coupled Schroedinger 
equations. The operator of the potential energy 
in Eq. (11) would then be 


w 
ir? rf 
Vio 

2m yw 


rn 


yw 

= 


which is positive definite if y>0, since both sum 
and product of its eigenvalues at any r are 
positive. If y=0, Vio is diagonal and not nega- 
tive. The dipole potential thus leads to no net 
attraction in the states ¢; and ¢: and no bound 
states of this type may be formed. 

The potential operator in Eq. (12) is 


ww w 
ht 
Vu=— 
2m 7 
r 


which is positive definite if y>v2. At y = v2, one 
eigenvalue of V,, is zero. The value y=v2 is 
obtained when F=1. However, F may also be 
zero, in which case y=0. Then V,; is diagonal 
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and has the single non-vanishing element 
w 


—— 


2m\r? 


which is everywhere positive only if w<0. If w>0, 
as in the case of neutron-positron interaction, 
there is a region of strong attraction between the 
roots of Vy, viz., =4$(2+V2)w. At r=w, the po- 
tential is —h?/2mw* and the range of negative 
potential is approximately v2w. Thus, the product 
of the potential and the square of the range is of 
the order of h?/m, as is required to form a bound 
system. It will be noted, also, that as r-0, Vis 
becomes positive and infinite so that there is no 
difficulty with the solutions at that point. 

The foregoing analysis shows that solutions 
of the wave equations may be expected to exist 
without altering the dipole potential near r=0. 
Of the possible solutions, the case of y =0 is most 
interesting because, in this case alone, there is a 
net attractive potential in the system. Under 
this condition, however, the equations for ¢2 
and ¢3 degenerate to those for plane waves, and 
the equation for ¢; contains a repulsive po- 
tential for all radii. Thus, a bound state must be 
characterized by ¢1=¢2=¢3=0. This condition 
is consistent with ¢,~0, according to Eq. (10), 
only if E=mc and 

1 w 
—-—— (13) 
dr r # 


Equation (13) can be integrated at once, up to 
a factor, 


1 
(14) 


That the solution, Eq. (14), represents a bound 
state is apparent from the fact ‘that the wave 
function, f=¢,/r, falls off as r~* with large r, so 
that most of the charge of the positron is 
concentrated within a radius of the order 
|w| =|»|e?/Mc*. The existence of such a state 
thus depends upon the dipole field persisting to 
radii small compared with the classical radius of 
the proton. This condition is not expected to be 
fulfilled in the case of the neutron because of 
the mesonic origin of its magnetic moment. But 
even if the range of the meson field were small 
compared with w there is serious doubt that 


effects connected with the time rate of change of 
the neutronic dipole moment would be 

in a state of such small radius. In particular, the 
inertial effects in the magnetic field might be 
expected to be large and, since these contribute 
a positive energy, the bound state would disap. 
pear into the continuum. 


2. SOLUTIONS FOR F=0 


General solutions, for which »?#0, are ob. 
tainable for the waves in which F=0. As stated 
above, only two wave functions are involved, 
¢1 and ¢,4, and the solutions are determined 
the equation of second order, Eq. (11) 


plus one of the relations, Eq. (10), say 


he 4 w i 


First, we transform Eq. (15) to a resemblance 
of Bessel’s equation by choosing a new dependent 
variable, y, 


"1 
At large r this approaches the equation for 
J=,(pr/h), the +}-order Bessel functions. If we 
write Eq. (17) with 1/r as the independent 
variable the form of the equation remains the 
same except that p?/h? is interchanged with —w’. 
The asymptotic solutions at small r are, there- 
fore, J+,(iw/r). We shall be particularly con- 
cerned with two forms of asymptotic solution, 
yi and ye, and their product Y=yyy2: 


(r/w)! sinh(w/r), 
y2>(r/w)! cosh(w/r), (18) 
Y—(r/2w) sinh(2w/r). 


The differential Eq. (17) can be put in more 
familiar form by transformation of independent 
variable from r to x, such that 


r= (19) 


y(x) cosas ) x) 0. 
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Equation (19) is the Mathieu equation in its 
conventional form, but with complex parameters. 
The general solution of the Mathieu equation 
has been found by F. Lindemann.‘ The method 
is based on the proof that there are two inde- 

dent solutions the product of which is an 
integral function of the variable z, which is 
related to x by 

2=cos*x. 


If Y(s) is the required integral function, the 
two solutions are given by 


where C is a constant such that 


An asymptotic expression for Y appears in 
Eq. (18) and this may be written as a function 


of z: 
for ’°’<hw/p, 


sinh4(i 

A(ipwz/h)i 
It is readily verified that if the product, Y, 
contained y;* or yz’, as given in Eq. (18), its 
asymptotic form would not be an integral func- 
tion of z. If we substitute Eq. (21) into Eq. (20) 
and integrate, we actually find 


(r/w) coth(—ipw/h)'} sinhw/r 
y2[(r/w) tanh(—ipw/h)*}! coshw/r 


so that the solution that is bounded and square 
integrable at r=0 must be of the form 


y=Y2—41 tanh(—ipw/h)}. (22) 


Since Y(z) is an integral function, the same 
asymptotic form of solution holds for all large z. 
If we apply the limiting process to Y, but such 


that 
rD>hw/p, 2—}i(pr*/hw), 


_ ‘F. Lindemann, Math. Ann. 22, 117 (1883). The method 
is described in Forsyth, Theory of Differential tions, 
Part III (Cambridge, 1902) p. 431 and also by P. Humbert, 
Fonctions de Lamé et Fonctions de Mathieu, Memorial des 
Sciences Mathématique (Gauthier-Villars, 1926). 


we get 


yi cos(pr/h), 
sin(pr/h), 


as 


Substituting these forms into Eq. (22) we get 
the asymptotic form at large r for the solution 
that is acceptable at r =0, viz. 


/h 
—tan(ipw/h)* tanh(—ipw/h)' cospr/h]}. (23) 


Since w corresponds to a length of the order 
e?/Mc*, the product pw/h is of the order unity 
only in the neighborhood of billions of electron 
volts. Such energies, however, are beyond those 
allowable in this treatment because the relativ- 
istic properties of the neutron have been neg- 
lected. In application, therefore, pw/h must be 
much smaller than unity and, disregarding the 
constant factor, the solution for y is essentially* 


yor? sinp(r—|w|)/h. (24) 


In brief, the electron waves manifest a phase-shift 
of —p|w|/h radians at large radii, if they form 
a state of total angular momentum zero, F=0, 
with the neutron. The same phase shift applies 
to both ¢; and ¢, components, as is evident from 
Eq. (16). The oS, and oP, waves at large r will 
be of the usual type, modified by this phase shift. 
This analysis holds for all finite ». When p is 
exactly zero the solution to the wave equations 
is the bound state presented in the preceding 
section. 


3. NATURE OF SOLUTIONS NEAR THE 
ORIGIN 


The behavior of the solutions of Eq. (10), or 
the equivalent Eqs. (11) and (12), near r=0 may 
be investigated by disregarding the zero power 
of r, ie., setting p?=0. The functions thus ob- 
tained will then be asymptotic forms of the 
complete solutions. Certain linear combinations 
of the asymptotic forms will be required, in 
general, in order to satisfy the boundary condi- 
tions at r=0. This study has been made for the 
states of lowest angular momentum, F=0 and 
F=1, which appear to be particularly tractable. 


*If w<0, the analysis holds with p<0 and the nega- 
tive square root in Eq. (22). 
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In fact it was shown in the preceding section 


that the case of F=0 is readily soluble. 


When F=0 we shall designate the asymptotic 
forms of ¢; and by and respectively. 
The functions ¢: and ¢3 are irrelevant. The 
equation for u:, as derived from Eq. (11) by 


setting =, and p’=0, 7 =0, is then 


w? 
=0. (25) 


This equation has the solutions, up to a factor, 
u=rexp—w/r and u,=rexpw/r. (26) 


One of these solutions is bounded, the other 
unbounded, at r=0 depending upon the sign of 


w. From Eq. (12) we get the equation for u, 
—4w/ +w?/ U4 =0, (27) 


which has the solutions 


(28) 


rT 


The first form is bounded and square integrable 


at r=0 if w>0, the second if w<0. The numerical 


factors with which u; and u, occur in the com- 


plete solutions are determinable from Eqs. (10). 
For P, waves we get directly from Eq. (10) 


uy’ +(w/r?—1/r)u, =0, 
ug’ + (1/7 =0, 


with acceptable solutions 


2 
(oP 
f 
0 r’ 


1 
u=-e'", w>0. 
r 


The first’ set holds if w<0. The second set is 


obviously the bound state described in Section 1. 


For S, waves, we replace mc in Eq. (10) by 
and set E=me’, ¢1=u1, d4= — uy, leading 
to 


uy’ +(w/r? 
ua’ 


with the solutions 


u,=re""", us=0, w <0, 
ug= —(hw/mc)-e-'", 
r 


The asymptotic forms of solution for F=1, 
i.e., y=V2 in Eq. (10) may be found by a series 
of transformations and quadratures with p? =0, 
The steps in this process will be omitted here 
and the results converted to those for a, b, f and 
g of Eq. (8). It is to be understood that these 
letters stand for solutions with p?=0 as used in 
this section. In particular, E is set equal to me 
in Eq. (8). There are four solutions in this case: 


(29) 
a= =U, =U, 


w w 
a= 
r 
(30) 
w w 
b=0, g=0, 
w? 2w* dr’ 
0 


4w? 
pon —(w? +wr) 
9r? 


w 13 
v2 9 9 


(31) 


SGP wy 


33 


QS Ras 


| 
| 
4 w dr’ 
+-u(1-—) f evlr’__ 
3 3r 0 r’ an 
4/w w . dr’ 
| f 
9 r? 0 r’ Ss 
5w h h 
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2w* 
Xe-wlr f 
r 
8w 
ac wrt +— 
9 9r 


(32) 


2w* 
f 
9r? r’ 


5 w\ h 
12 r 


If w>0 the solution of Eq. (31) is unacceptable. 
Equation (32) defines a D-wave that is unique 
except for an undetermined lower limit, ro, in 
the integrals. This lower limit must be finite, 
however, and admixes the solution (30) in terms 
of the order w*. Both solutions (29) and (30) are 
acceptable S-waves. It is evident, consequently, 
that the solutions are not completely defined 
and that additional boundary conditions should 
be imposed at r=0. Such boundary conditions 
might be stated by defining the central dipole 
more closely as the limiting case of approaching 
magnetic charges or as a shrinking amperian 
current. This was first pointed out by Bloch® as 
having a major influence on the total scattering 
cross section. Its effect thus evidently appears 
in first order (in w) in states of higher angular 
momenta, but does not appear at all for F=0. 
In the above solutions for F=1 the uncertainty 
appears only in second, and higher, power of w. 


To first order, Eq. (29) and Eq. (30) are identical 


and the two independent solutions may be 
written 


’ g—0, (33) 


*F. Bloch, Phys. Rev. 51, 994 (1937). 


D: oer, 


If w<0, Eq. (30) is unacceptable, but Eq. 
(29) becomes a unique S-wave. Both (31) and 
(32) provide D-wave solutions, if ro=0, and 
they are identical to first order in w. The approxi- 
mate solutions in this case are then those in 
(33) and (34), as before. 

In addition to the ;S,; and ,D3,;2 waves there 
are the 1P, and :Ps/2. The equations for the 
latter may be obtained essentially from Eq. (8), 
or their transformations Eq. (10), by replacing 
mc by —mc*. The second-order equations (11), 
which determine } and g, are unchanged in the 
transformation. The solutions for 6 and g pre- 
sented in (31) and (32) can thus be combined to 
give solutions for ¢; and ¢2. These solutions are 


(35) 
gi=—(2r?+wr)v2, —(2r?+2rw+w*). 


However, Eq. (11) must have four sets of solu- 
tions. The remaining two sets can be derived 
from (35) by noting that Eq. (11) is invariant to 
the simultaneous substitutions w— —w, ¢:— — ¢1. 
Hence two new solutions are obtained from (35) 
—V2 rel", (36) 
o2= — (27° —2rw+w*). 
Again, the wave functions are unambiguous 
only to first order. The result may be written: 


v2 
br — jw, 


“(3 wv2 h (37) 
2 mc 4r mc 

v2 
1P 3/2: — g—r+ Fu, a—0, f-0. (38) 


4. SOLUTIONS TO FIRST ORDER IN w FOR F=1 


Equation (8) may be expressed as two sets of 
coupled, second-order differential equations. 
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When F=1 these equations are 


a’ w 
a” +2—+—a+v2—f -—(a+v2f) =0, 
r 


et =0, 


bo fp? 4h 
b’+2—+( —-— (39) 
r 3r8 
0, 
3r4 
+2 
=0. 
3r8 


If w=0, these equations become four independ- 
ent differential equations closely related to 
Bessel’s equations of order, 4, 3/2, and 5/2. The 
solutions of Eq. (39) that are bounded at r=0 
are then, to zero order in w 


ao= 
bo = 
go= (rh/2pr) 
fo= 


(40) 


where ao, Bo, Yo, and 4» are constant factors. 
Expanding the solutions of Eq. (39) in an 
ascending power series in w, the functions Eq. 
(40) appear as the leading terms. The functions 
that are multiplied by the first power of w may 
be found by substituting the functions of Eq. 
(40) in those terms of Eq. (39) that appear to 
the first power, neglecting the terms in w* and 
solving the resulting, inhomogeneous differential 
equations. The integrals will consist of two parts, 
a particular integral and a complementary func- 


tion ; the latter is chosen to be orthogonal to the 
solutions Eq. (40) at large r. Let ay, 6;, 1, and 
5, be the factors with which the complementary 
functions in a, }, g, and f are added, » Tespectively, 
The solutions to first order are then computed 
to be 


4 
eu (=) [aos 
v2 
ré h 


— 
b= (=) [+07 3/2(€) +81J_3/2(€) 


4Bot+v270 w 


g= (=) [rw +1J_3/2(€) (41) 
—4y0 w 
i(é 
J 
f= (=) oJ 5/2(€) 


w 
3/2(&) —| 450+ 
r? 


w 
(vo—4V260) 
3ré 


Reversal of the sign of mc*, as is done for half of 
the solutions changes only the form of f as 
evidenced in the explicit appearance of mc*. For 
the rest, the numerical coefficients in Eq. (41) 
may be determined to fit onto the asymptotic 
forms near r=0, as presented in Section 3. 
Making use of the asymptotic forms of the Bessel 
functions® for large pr/h, the induced phase 
shifts are brought into evidence. 

The asymptotic forms, at large r for the 
solutions Eq.: (41) that fit onto the accepted 


6 Jahnke Emde (B. G. Teubner, Leipzig, 1933), second 
p. 203, for example. 
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behavior near r=0 are as follows: 


r 
w 
b—- ) cosp(r-—) 
r\E+mc 3 


v2 pw 
cospr/h, 


v2 E-me\3 
) sinpr/h, 


12 hr \E+mc? 


3/2 


(42) 


v2 pw 
a—>— —— COS h, 
12 hr 


— sinpr/h, 

and the related expressions for ,P; and ,P3;2 
which are obtained by exchanging @ with —), 
f with g, and sin for cos in the expressions for 
»Ss and ,D3,2, respectively. The results of the 
above tabulations may be summarized as follows: 
electron states forming a state of total angular 
momentum hf, with the neutron, show a phase 
shift of 4pw/h radians if 7=3/2 and —4}pw/h 
radians if j=}. In addition, each wave is at- 
tended by the factor +(v2pw/12h) times the 
unbounded wave belonging to the alternate state 
of the same j. This is to be compared with the 
phase shift —p|w|/h radians for oP, and oS, 
states. 


5. CONCLUSION 


The foregoing results for the phase shifts 
indicate that the scattering of an electron by a 
magnetic dipole is characterized by a cross 
section of the order of rw*. This has been pointed 
out already by Massey.’ For the electron-neutron 


7H. S. W. Massey, Proc. Roy. Soc. A127, 666 (1930). 
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interaction, therefore, the cross section is of the 
order of 10-*° cm*. The theory predicts, however, 
that under certain circumstances the amplitudes 
of waves scattered from many. neutrons in a 
nucleus should reinforce even though the spins 
of the neutrons are distributed statistically. This 
appears to be the case for S-waves, for example 
(the same argument applies to Py-waves), since 
the phase shift in 9Sj-states is —p|w|/h and 
that in ,S, states is —pw/3h. The weighted 
average phase shift is then —}pw/h if w>0 and 
zero if w<0. Thus, if the magnetic moment of a 
neutron is localized sufficiently that each neutron 
in a nucleus of N-neutrons acts as an individual 
scattering center, and if the wave-length of the 
incident particle is long compared with nuclear 
diameters, the cross section for positron scatter- 
ing (w>0) in S-waves should be of the order of 
N*w*, whereas that for electron scattering in S- 
waves would be of the order w* and would 
depend upon the resultant magnetic moment of 
the nucleus. 

Similar conclusions to those found for the 
neutron-positron interaction probably can be 
drawn for the electron-proton system. In partic- 
ular, the strong attractive potential in the oP,;- 
state may have some effect upon the term values 
of the hydrogen atom that would illuminate the 
consequences of some of the assumptions that 
have been made. In general, the motion of an 
electron in a combined electric monopole and 
magnetic dipole field should be of more funda- 
mental interest than that in the magnetic field 


alone. Such solutions, especially if generalized to 


include unit spin for one of the particles, would 
apply to the interaction of the ‘“elementary”’ 
particles for which electric charge and magnetic 
moment are necessarily associated. On the other 
hand, the chief interest in such solutions prob- 
ably would center on the possibility of bound 
states and the theory of such states is rather 
clouded by the uncertainties due to spin inertia 
and effects of retardation. 

The work reported in this paper was done 
during the author’s association with The George 
Washington University and The Department of 
Terrestrial Magnetism, Carnegie Institution of 
Washington, and it is a pleasure to acknowledge 
the support of these institutions. 
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Note on the Scattering of Radiation in an 
Inhomogeneous Medium{ 
C. L. PEKERIS 
Columbia University,* New York, New York 
January 23, 1947 
E consider the scattering of a plane sound wave 
when passing through a medium in which the sound 
velocity ¢ fluctuates slightly and irregularly around the 
mean value é. The corresponding problem of the scattering 
of electromagnetic radiation in an optically inhomogeneous 
medium can be treated along similar lines. We write 
c(x, 2) =t+Ac(x, 3), (1) 
and compute the intensity of sound scattered by a volume 
V of the medium, which is large enough to contain several 
oscillations in sign of Ac. We choose a coordinate system 
centered at V and assume an incident plane wave traveling 
in the negative z direction. If g denotes the sound potential 
then we have 


i= 
where 
du=dx'dy'ds’, (3) 
and c has been written for ¢. From an analysis of the 
scattering by a sphere of radius a of slightly different 
sound velocity it is found that expression (2) is valid if 
_ (ak)(Ac/c)<“«1. We should, therefore, expect the perturba- 
tion method used in deriving (2) to be valid if (a’k)(Ac/c)<1, 
where a’ denotes the dominant “‘wave-length” of Ac. 
If the distance r of the point of observation is large in 
‘comparison with the wave-length, we have 


where we have taken the real part of (2). The energy E 
scattered in an elementary solid angle dQ is given by 


(5 


NUMBER 4 


We now average over time, ont writing Eo= tock, obtai 


E=E (20) aw 


‘Ac 
ov 
xf 9 +8, 
xB 
(6 
with 
B=ay"—y’, y=s"—2'. 
system, and the limits of a, say, are —(x:+23) and 
+ (x:+<2) if the limits of x are —x, and +22. In particular 
we have 
(1 +5)a 
+sino sin sin(y— go), (7) 


where 60, yo refer to the direction of r, and @, ¢ are the 
spherical coordinates in the (a, 8, y) system. We now define 
the autocorrelation R(a, B, y)* of the values of (Ac/c) at two 
points inside V whose x’, y’, 2' coordinates differ, respectively, 
by a, B, 7 by: 


= Rla Br). 
Substituting in (6) we get 


6) 
In the case of spherical isotropy we have 
R(a, 8, (10) 
and using the relation 
+sin@ cos(y— yo) de 


we obtain 


For a Markoffian autocorrelation function* 
R(p) (13) 


where a represents a measure of the “‘size”’ of the scattering 
elements, and Eq. (12) yields 


It should be noted that 69 is measured from the negative s 


a=x -x, 
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axis. If the angle is measured from the forward direction of 
the incident wave, cos(@0/2) in (11) and (14) should be 
changed to sin(@0/2). 

In the interpretation of light scattering in solutions, it is 
customary to assume a model for the scattering elements, 
say a sphere of radius a, and apply the theory of scattering 
by a cloud of randomly distributed spheres.* Experiments 
indicate, however, that a has a considerable statistical 
scatter. The model of randomly distributed spheres of 
variable radii is perhaps not simpler to visualize than the 
autocorrelation function R itself, which is the only. sta- 
tistical property of the scatterers’ population determining 
the angular distribution of scattered light, as shown in 
Eqs. (9) and (12). 

This note is taken from a study of scattering problems 
made by the writer in May, 1942 for O.S.R.D., of which 
only the results were published. The manuscript was made 
available to members of Division 6 and was loaned to 
several members of Division 17 of O.S.R.D. Among the 
problems treated in this study, either exactly or approxi- 
mately, are the scattering of sound by (1) a vortex sheet, 
(2) a cylindrical vortex, (3) a spherical vortex, (4) a 
convective cloud in which turbulence is isotropic, and (5) a 
surface. 


tion assisted by the Ernest Kempton Adams Fund for 
of Columbia University, New York, New York. 
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The Half-Life of Be” 


D. J. HuGues, C. EGGLer, anD C. M. HuDDLESTON 
Argonne National Laboratory, P. O. Box 5207, Chicago 80, Illinois 
January 9, 1947* 


HE long-lived nucleus Be” has been reported as 
formed only in the Be*(d, p) reaction.+* The half-life 
could only be estimated to be in the range 10* to 10’ years'~* 
because the yield of the reaction producing the Be” is not 
known. Bretscher‘ tried to produce Be” from the Be*(n, ) 
reaction but found no activity. Levinger and Meiners? also 
found no Be” activity in Be irradiated in the Clinton pile. 
Because the extremely long half-life of Be" is of interest 
theoretically, it was decided to attempt its measurement 
by isolation of the activity produced by the (, y) reaction 
in two BeO samples which had been irradiated for different 
periods of time in the Hanford pile (kindly supplied by 
O. C. Simpson of the Argonne Laboratory). Although the 
BeO was of high purity, it had a gross counting rate of 
about 5000 counts per min. per mg of Be. As the expected 
Be” activity would be only a few counts per min. per mg, 
it was necessary to purify the samples carefully. The 
purification was carried out by group separation and the 
addition of suspected impurities of the third and fourth 
groups as carriers, the final specific separation for Be 
being the basic acetate process. The main impurities were 
found to be carried with titanium in excess NaOH and 
with zinc in the H2S-NH,CH;COO process. 
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The final samples, which did not decrease in activity with 
further purification, had activities of 3.3 counts per min. 
per mg of Be for one sample and 6.7 for another sample 
that had been irradiated twice as long. Absorption curves 
of the activities gave an energy of 0.58+0.03 Mev which is 
in close agreement with the findings of Levinger and 
Meiners? and McMillan and Ruben! and the mass differ- 
ence of Be” and B” (0.58 Mev’). 

In order to obtain the half-life from the mvt and observed 
activity it is necessary to know the cross section of the 
(n, y) reaction. The cross section for Be*(n, y) was meas- 
ured by comparing the effect on the reactivity of the 
Argonne graphite pile of several kilograms of BeO with that 
of graphite. The value obtained was 0.01 barn (barn 
=10-** cm*) in agreement with the measurements of 
S. K. Allison,* of A. Langsdorf,* and of R. Nobles.* The 
half-life of Be” calculated from the activity, mvt, and the 
cross section is almost the same for the two samples, 
the wae value being 

Ty =2.9X 108 years. 


The amount of activity to be expected in the Hanford 
graphite from the C¥(n, a)Be” (fast neutron) reaction was 
then calculated, using the measured half-life and an (n, a) 
cross section estimated by comparison with measured pile 
cross sections for other endothermic reactions.’ A small 
amount of graphite which had been in the Hanford pile 
for about a year was ignited, several mg of Be carrier 
added, and the Be purification carried out. An activity was 
obtained of magnitude within a factor of two of that pre- 
dicted, which did not decay, and had the correct beta- 
energy for Be”. It seems certain then that the activity was 
Be” produced by the C(m, a) reaction and the finding 
constitutes additional support for the correctness of the 
half-life determination. 

Because of the importance of obtaining Be” in large 
amounts and of higher specific activity than can be done 
by the (m, y) reaction, arrangements are now being made 
to irradiate boron in the Hanford pile for the production 
by the (m, p) reaction and subsequent chemical separation 
of the Be”. 

Lavine and E. Meine Pt Phys. 

: Seaborg, Rev. Mod. Phys. i 16, 1 (i944). 

‘E. iauiee Nature 146, 94 (1940). 

* Noyes and Bray, A System of itative Analysis for the Rare 
Elements (The Macmillan Company, New York, 1943). 


Hughes, Manhattan | work (1943-46). 
N. Goldstein, and W. D. B. Spatz (unpublished). 


Rev. 70, 123 (1946). 
to be published. 


Photoelectric Emission from a Rough Surface* 


HuGH BRADNER 
Radiation Laboratory, University of California, Berkeley, California 
January 10, 1947 


XPERIMENTS on the normal energy distribution of 
photoelectrons from alkali metals ordinarily make use 
of evaporated metals which have a matte surface rather 
than the smooth surface considered in the Fowler-Du- 
Bridge’ theory. The following is a modification of the 


| 
! 
ain 
| 
) 
| 
7) 
y, 
2C,. L. Pekeris, Phil. Mag. , 541 (1942). 
+P. Debye, J. App. Phys. 15, 338 (1944). " 
| 
) 
) | 
) 
| 
) 
) 


270 LETTERS TO 


At 0°K, the electron velocity distribution is uniform 
within a sphere of radius (240/m)+ in velocity space, where 


#o= maximum energy of unexcited electrons, and m=elec- 


tron mass. Assume that the energy of the photon adds to 
the electron energy in such a way that it increases only the 
component of electron velocity normal to the surface 
element. And assume further an excitation probability a 
proportional to the normal component of the electron 
velocity. Considering the transformations of the sphere by 
the excitation energy hy, the surface potential barrier Wa, 
and the retarding potential V, applied at an angle 6 to a 
surface element, one obtains an equation that is readily 
soluble when sin*é. The current from 
an element is 
i= (2rema/h*) cosé- [uothy — W.- VeF, 

which is just DuBridge’s expression multiplied by cosé. 

The photoelectric current from a rough surface at 0°K 
thus shows the same dependence on retarding potential 
and on energy as that from a smooth surface, providing 
the surface is not too rough, and the retarding potentials 
are not too low. Assuming a maximum angle of 45° between 
surface element normals and the normal to the cathode, 
the formula is valid for i=(2+wema/h*)(V.)*. 

The equation for V.<(ue+hv— Wa) sin*d does not re- 
duce to a simple form. From the equations it does follow, 
though, that the current will be less than DuBridge’s 
expression, indicating that photoelectric saturation current 
curves for rough surfaces will be rounded even at 0°K. 

The writer has examined the data of Overhage* on normal 
energy distribution to see whether the above arguments 
could explain any of the discrepancy between theory and 
his curves. They do not. 

The writer wishes to express his indebtedness to Dr. W. 
‘V. Houston who suggested looking into this problem, and 
also suggested the mathematical treatment used. 

Overhage, Phys Rev. $2, 1099 41937). 


Neutron Induced Activities in Lutecium, 
Ytterbium, and Dysprosium 


Mark G. INGHRAM, RICHARD J. HAYDEN, AND Davin C. Hess, Jr. 
Argonne National Laboratory, University of Chicago, Chicago, Illinois 
January 17, 1947 
ITRIC acid solutions of Lu,O;, Yb,0;, and Dy,O; 
were irradiated with neutrons in the Argonne heavy 


water pile. Decay curves of the irradiated lutecium sample 


showed 3.4-hour and 6.6-day activities as reported by 
Flammersfeld and Mattauch.' Decay curves of the ytter- 
bium sample showed a 2.7-hour activity (probably the 
same activity reported by Marsh and Sugden? as 3.5 hours) 
and also a 102-hour activity not previously reported. No 
evidence of the 41-hour activity reported by Pool and Quill* 
was found. Decay curves of the dysprosium sample showed 
the 2.5-hour activity reported by Marsh and Sugden.’ 
Assignment of mass to three of these active isotopes was 
made as follows. A portion of the sample containing the 
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activity in question was placed on the filament source of 
a mass spectrograph. Operation of the spectrograph sepa. 
rated the isotopes and deposited them on a Photographic 
plate. The position of the active isotope on this plate was 
determined by placing a second photographic plate face 


to face with the original. The disintegration of the act; 


isotopes on the first plate produced developable images on 
both plates. Comparison of the position of this active 
deposit with the normal mass spectrum on the original 
plate determined the mass of the active isotope. This 
method of mass assignment gives results without the 
necessity of assuming the nuclear reaction by which the 
isotope is formed and is therefore frequently valuable in 
establishing the nuclear reaction. Using these techniques 
the following masses were assigned: 6.6-day lutecium, 177. 
102-hour ytterbium, 175; and 2.5-hour dysprosium, 165, _ 
1 Flammersfeld and Mattauch, Naturwiss. 31, 66 (1943), 


2 Marsh and S Nature 136, 102 (1935). 
3M. L. Pool L. L. Quill, Phys. Rev, 53, 437 (1938), 


Interpretation of Anomalous Larmor Frequencies 


in Ferromagnetic Resonance Experiment 


CHARLES KITTEL* 
of Electronics, Massachusetts Institute 
, Cambridge, Massachusetts 
January 10, 1947 


ECENTLY J. H. E. Griffiths' reported an important 
new ferromagnetic resonance experiment at micro- 
wave frequency. The experiment is essentially the ferro- 
magnetic analog of the Purcell-Torrey-Pound nuclear 
resonance experiment. Griffiths found however the unusual 
result that the resonance frequencies he observed were 
greater than the calculated Larmor frequencies for electron 
spin by factors of about two to six. He attempted unsuc- 
cessfully to explain the anomaly by the introduction of the 
Lorentz cavity force, a step which is definitely not justified, 
A theory of the resonance effect is given below which 
leads to values of the resonance frequency in close agree- 
ment with the experimental determination. It is shown that 
it is important to consider the dynamic coupling caused by 
the demagnetization field normal to the surface of the 
specimen. The result is that the appropriate Larmor fre- 
quency should be calculated as for a field strength (BH) 
rather than for H. 
A ferromagnetic specimen whose surface is the plane 
y=0 is subjected to a strong d.c. field H, and a weak 
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' microwave field H,. The magnetization M and the angular 


momentum density J are related by M=yJ, where 
y is the gyromagnetic ratio. The equation of motion 
8J/at = XH] may be written 

8M/at=y[MXH], (1) 


where the components of H are (Hz, —4xM,, H,). Here 
—4nrM, is written for H, on the basis of the divergence 
relation B, = H,+42M,=0. The Lorentz and Weiss fields 
are omitted as they are always parallel to M, and hence 
their vector product with M is identically zero. It is as 
sumed that saturation obtains, and in fact that the dc. 
field is sufficiently strong to outweigh the magnetic 
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TaBLE I. Comparison of observed and calculated resonance frequencies. 


(rad /sec.) 
d.c Calculated Larmor Experi- 
field frequencies mental 
for H, for (B.H:)§ frequencies 
Fe 2800 5.0 14.5 15.4 
500 0.9 5.8 5.9 
Co 510 0.9 5.3 5.9 
Ni 5000 8.8 13.5 15.4 
3800 6.7 10.9 13.2 
1030 1.8 4.9 5.9 


anisotropy forces so that in the static case the spin system 
always points in the direction of the resultant magnetic 
field. In this limit the whole specimen is supposed to behave 
as a single domain, so that there are no complications from 
the domain boundaries. 

We have, neglecting products of small quantities, 


0M, =7B.M,; (2) 

(3) 

aM,/at=0; (4) 

which gives —w*M,=~7°B,(M.H,—M-H,). The suscepti- 
bility xe = M./H, is given by 

xs = (S) 


where the frequency for resonance is found to be wo 
=~+(B,.H,)' and is the Larmor angular frequency in the 
fictitious field (B,H,)*. The static susceptibility xo is equal 
to M,/H, and varies with the biasing field; here M, may be 
taken as the saturation magnetization. 

In Table I Griffiths’ experimental values of wo are com- 
pared with the calculated Larmor frequencies for fields H, 
and (B,H,)'. The gyromagnetic ratio is taken as the 
electron spin value. The induction B, is calculated using 
1700, 1400, and 500 as the saturation magnetization for Fe, 
Co, and Ni, respectively. The theory given here is seen to 
account quite well for the experimental results, considering 
that damping and anisotropy forces have been neglected. 

The effect of relaxation is introduced by adding a term 


—h(M—xoH) to the right side of Eq. (1). We find 


wo? yo +jwr (6) 


where yo= B,/H,. From this it is seen that the demagnetizing 
field acts to increase the apparent damping. 

The interpretation of a complex susceptibility in terms of 
the result of an electrical measurement is somewhat 
involved, but may be accomplished by the reasoning 
previously given.? The apparent permeability for a re- 
sistive measurement (such as a “Q" measurement) is 
ur=|uz|(1—sing), where us is the permeability derived 
from xz, and ¢ is the phase angle of us; ¢ is always between 
Oand —r. 

* John Simon Memorial Foundation Fellow. 


Guggenheim 
1J. H. E. Griffiths, Nature 158, 670 (1946). 
*C. Kittel, Phys. Rev. 70, 281 (1946). 
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Biased Betatron in Operation 


W. F. WEsTENDORP 
Research Laboratory, General Electric Company, Schenectady, New York 
January 22, 1947 


BIASED betatron'~ with closed central core has been 

successfully operated. The guiding magnetic field 
consists of a direct-current field component with a sinus- 
oidal component superimposed, while the closed central 
core carries only a sinusoidal component of flux main- 
tained at the proper magnitude by the current in a bucking 
coil located in grooves in the polefaces. Figure 1 shows a 
schematic cross section of the machine and Fig. 2 is a 
diagram of the principal electrical components of the 
energizing circuit.! The turn ratio of main coils to bucking 
coil is chosen so that the proper flux ratio and thereby the 
betatron 2:1 rule is fulfilled. Although radically different 
from previous betatrons in which the 2:1 rule was estab- 
lished by means of air gaps, this principle of locating the 
electron orbit offered no difficulties. It has the advantage 
of a reduction of the size of the capacitor bank since less 
magnetic energy storage takes place in a machine with a 
closed central core. The application of a direct-current 
bias to the main coils to produce an increased guiding 
field and to the bucking coil to keep the constant com- 
ponent out of the central flux, proved to be not critical and 
was tried out successfully for all ratios between main coil 
direct-current and crest alternating-current between 0 and 
0.866. The corresponding ratio of current components in 
the bucking coil was always smaller for maximum x-ray 
output, but not critical. 

No compensating or phase correcting circuits of any 
kind were used. Even though the a.c. components in main 
and bucking coils are not exactly 180 degrees out of phase,* 
theory*® predicts that the shift of the orbit produced is 
very small during most of the accelerating interval and 
this was checked experimentally by the use of an orbit 
tilt circuit by means of which the electrons are made to 
hit a target located above the orbit. 


CAPACITOR: 
BYPASS. 
REACTOR 


GROOVE 
CONS 
0.C. CONTROL 
= RESISTANCE 
aii. 
0.c. 
GENERATOR 


Fic. 1, Biased betatron with closed central core. 
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Fic. 2. Circuit diagram of biased’betatron. 
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By means of a “megavolt meter,‘ a direct reading of 
the electron energy and of the corresponding x-rays pro- 
duced could be obtained in all cases. Thereby an instan- 
taneous change from 11 Mev without d.c. bias to 20.5 Mev 
with d.c. bias could be produced and read by merely turn- 
ing on the direct current and rephasing the electron injec- 
tion pulse. In both cases, electron injection voltages as 
low as 8000 volts could be used indicating a high degree of 
uniformity of the magnetic guiding field with azimuth. 
With oil-cooled coils, the machine will produce 50-million 
volt x-rays. 

iw. F. Kerat Phi, Rev. 68 Phys. 16, 285 (1945). 


E. Amaldi and B. Ferretti, Rev. beg oe ag , 389 (1946). 
We F. Westendorp, Rev. Sci. Inst. 17, 218 lidiey. 


iuitiiinis Cross Section of Boron for Neutrons 
of Energies from 0.01 ev to 1000 ev* 


R. B. Sutton,’ B. D. McDantez,? E. E. ANDERSON,® 
AND L. S. LAVATELLI* 


Univer. Ala Scientific Laboratory, 
ae Santa Fe, ‘New Mexico 


January 23, 1947 


HE slow neutron capture cross section of boron has 
been measured for neutron energies from 0.01 ev to 
1000 ev. The cross section was obtained by determining 


' the transmission of boron-containing samples. Neutron 


energies were determined by using the method of modula- 
tion of a cyclotron beam to obtain the time of flight of 
the neutrons over a 7.6-meter path between the neutron 
moderator and a BF; neutron detector. This method is 
similar to that described by Baker and Bacher® and by 
Bacher, Baker, and McDaniel.* The equipment used was 
constructed in the Physics Department at Cornell Uni- 
versity by B. D. McDaniel and C. P. Baker; a description 
will appear in the literature shortly. Neutron pulses in 
twelve time-of-flight intervals could be recorded simul- 
taneously. 

In order to cover the range of neutron energies four 
samples, each of different boron content, were used. The 
one containing least boron was a BF; gas cell; the others 
were B,C. Of the three B,C samples the thickest and 


thinnest consisted of cells containing B«C powder; the . 


third was composed of a mixture of B,C and lead borate 
which was hot-pressed to form a compact disk. 

The curve of total cross section or vs. time-of-flight 7 for 
each absorber was found to vary in the following manner 


or=ar+8, 


where a and £ are constants. It was assumed that the 
constant 8 corresponded to the scattering cross section of 
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the various elements present, and the term ar 

to the absorption cross section of boron. The sca 
cross section, 8, was obtained by taking the pie 
value of ¢ at r=0. The scattering cross section for boron 
was found by correcting 8 for the cross sections of the other 
elements present. Values used were for fluorine, 4.1) 10-« 
for carbon, 4.85 X 10-*4 cm*, for oxygen, 4.2X 10-* om? 
and for lead 9.6X 10~** cm*. The value obtained for the 
boron scattering cross section was 4.2 X 10-*4 cm?. The total 
contribution of the other elements in the B,C samples was 
2.2 10-** cm? per B atom. The capture cross section of 
boron was obtained from 


a=o7r—B=ar 


by assuming that all the capture cross section was repre- 
sented by ar and none was included in 8. In Fig. 1 the 
values of the product of the capture cross section by 
neutron velocity are plotted as a function of neutron 
energy. The instrumental resolution function is triangular, 
The resolution used depended on the energy region under 
investigation ; in Fig. 1 the base of the triangle is equal in 
length to the separation of the points. 


The value of ov determined from the weighted average | 


of the points on the curve is 1.61+.02< 10-** cm? meter/g 
sec. It can be seen from the curve that systematic devia- 
tions occur between different absorbers. This is probably 
owing to errors in the determination of the boron content 
of the absorbers. The value of ov determined from the BF; 
gas alone is 1.58 10~** cm? meter/y sec. which may be a 
better value than the one given above, since the boron 
content of the BF; was probably better known. Both these 
values are lower than the figure of 118 10-**E-+ cm? (or 
1.63 10-** cm? meter/y sec.) recently given by Rain- 
water and Havens,’ but are higher than the value of 
10-*4 cm? meter/y sec. obtained by Bacher, Baker 
and McDaniel.* All agree within the uncertainties given. 
The figure also indicates that av is essentially constant 
(to +10 percent) up to 1000 ev, if the assumption is 
correct that 8 contains only components of scattering cross 
section. 

* This work was carried out under contract between the University 
of | Rey and the Manhattan District, Corps of Engineers, W: 

* Now at Milwaukee-Downer Co! » Milwaukee, Wisconsin. 

4 Now at Harvard University, Cam ridge, Massachusetts. 

5C. P. Baker and F. 


. Baker, and D. McDaniel, Phys. Rev. 69, 


443 
7 RS ener and W. W. Havens, Jr., Phys. Rev. 70, 136 (1946). 


Remarks on Dr. Bhabha’s Note “‘On the Expan- 
sibility of Solutions in Powers of the 
. Interaction Constants” 


Niets ARLEY 


Institute of Theoretical Physics, Copenhagen, and Palmer 
Laboratory, Princeton, New Jersey 


January 13, 1947 
N a recent note Bhabha! has pointed out that the 
difference between the physical and the non-physical 
solutions in the classical electron theory proposed by 
Dirac? may be characterized by the fact that the former is 
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continuous in the interaction constant e at the point e=0, 
whereas the latter is not, and that consequently the former 
may be expanded into a power series in ¢, but the latter 

may not. Since the usual perturbation methods of solving 
problems in any quantum field theory which is the quan- 
tization of a classical field theory just consists in expanding 
into a power series of the relevant interaction constant, A, 
Bhabha, therefore, lays down, not as a property to be 
proved, but as a general postulate that we have to demand 
all physical solutions to vary continuously for \—>0 in 
order that they be expansible into a power series of \ 
from the point \=0. 

Firstly, continuity is, however, not a sufficient condition 
for expansibility into a power series. For example the 
function 7 Ind is both continuous and has a continuous 
derivative at the point A=0, but is, nevertheless, not 
expansible into a power series of \ from \=0, this point 
being an essential singularity since the function has here 
an infinite branch point in the complex \-plane. The cor- 
rect condition for expansibility is the much more severe 
condition that the solutions be analytic in A at \=0. 

Secondly, in the quantum perturbation method used for 
calculating transition probabilities, the equations to be 
solved are infinite systems of simultaneous linear differ- 
ential equations. It is true that the usual perturbation 
method consists in expanding into a power series of 4, 
thus tacitly assuming the solutions to be analytic in A. This 
assumption need, however, by no means be true even in 
the simplest case in which all the coefficients are constant 
throughout the complex plane. In a previous paper*® we 
have given a general mathematical theory of such dif- 
ferential equations and have shown that all forms of 
‘pathologies’ may occur, such as, for instance, equations 
having no solutions at all, equations having infinitely 
many solutions for the same initial condition, equations 
with solutions which have a continuous derivative which is 
not differentiable for any value of the variable, and so on. 
In particular, there exist equations with constant coeffi- 


’ cients which have a uniquely determined solution which 


may be found explicitly, but for which the usual perturba- 
tion procedure gives a divergent result even in the second 
approximation. (In this example the explicit expression 
shows that the solution contains just a term of the form 
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Ind mentioned above.‘ Now it is in fact possible to 
develop other methods of perturbation, such as expansion 
in series which are only partly power series in A, the 
coefficients containing also 4, which at any rate converge 
in this specific example. We are investigating these ex- 
pansions more closely and the results will be published 
elsewhere. 

Notwithstanding the fact that some of the divergence 
difficulties of the present quantum field theories, especially 
those connected with the self-energy problem, are certainly 
of a deeper nature which will require a quite new theory, 
presumably consisting in introducing Heisenberg’s funda- 
mental length /» in a suitable way,' it is, however, likely 
that some of them, especially those connected with transi- 
tion probability problems, may be removed within the 
present mathematical formalism by better perturbation 
methods, as is indicated by the damping theory of Heitler 
and his co-workers, for instance. To insist in principle, as 
proposed by Bhabha, on the present special form of 
perturbation method is, therefore, neither justified nor 
desirable, and to do so may even be a hindrance for further 
progress. 


Rev. 70, 759 (1946). 
Proc. Roy. Soc. A167, 148 (1938 
Arley Borchsentus, Acta Math. 76, 261 (1 ). 
w reference 3, p. 300. 
Heisenberg. Ann. d. Physik 32, 20 (1938). 


Erratum: Expanding Universe and the 
Origin of Elements 
(Phys. Rev. 70, 572-573 (1946)] 


G. GamMow 
The George Washington University, Washington, D. C. 


HE value of the space curvature entering into the 
formulae (1) and (3) must be considered to be 
expressed not in centimeters but in the units of the selected 
length /. Thus estimating the radius of curvature from 
the expression (3) we get the value of 1.7 10-"(—1)*x/ 
=1.7X10~(—1)! cm, or about two billion imaginary light 
years, instead of 0.2 imaginary light year as given in the 
Letter to the Editor. This does not change, however, the 
further arguments. 
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MINUTES OF THE MEETING AT Los ANGELES, CALIFORNIA, JANUARY 3-4, 1947 


HE 276th meeting of the American Physical 
Society was held on Friday and Saturday, 
January third and fourth, 1947, at Los Angeles 
in Harris Hall of the University of Southern 


California. This was our first meeting at that. 


University, which furnished excellent accommo- 
dations and management for our sessions. Our 
thanks are due particularly to R. E. Vollrath and 
G. L. Weissler. The attendance at the most 
populous session amounted to some 200. The four 
invited papers were admirable in content and 


presentation. The vice president of the Society 
was prevented by the grounding of his plane from 
presiding at the session of invited papers, and 
was replaced as Chairman by C. S. Van Atta. 
The other Chairmen were R. E. Vollrath and 
J. Kaplan. 


J. KaPLan 

Local Secretary for the Pacific Coast 
University of California 

Los Angeles, California 


Invited Papers 


Interferometry. W. E. Wittiams, Pasadena. 


Design of a Small Nuclear Reactor. Ropert Curisty, California Institute of Technology. 
Recent Progress in Applied Absorption Spectroscopy. R. R. Brattain, Shell Development Company. 
Recent Progress in Applied Mass Spectrometry. D. P. STEVENSON, Shell Development Company. 


Contributed Papers 


Al. The Atmospheric Bands of Oxygen. JosePH KaPLan, 
University of California, Los Angeles—A remarkable new 
afterglow has been discovered in which the well-known 
atmospheric bands of O; are emitted. This is the first time 
that these bands have been observed in emission. The 
bands of the system which were most readily photographed 
are the 0,0 band at A7594 and the 0,1 band at approxi- 
mately 48650. The latter band has not been previously re- 
ported in either absorption or emission. Using small dis- 


- persion and Eastman IN plates, it was possible to obtain 


a good picture of the 7594 band in a one-hour exposure of 
the afterglow. Several other bands of the system and several 
apparently new bands of oxygen were obtained by using 
longer exposures. The new bands have not yet been identi- 
fied. The exposure times were short compared with those 
needed to obtain comparable spectra of the auroral after- 
glow in nitrogen. The source of this new afterglow was a 
tube in which the auroral afterglow in nitrogen had been 
excited, and in which most of the nitrogen had been cleaned 
up. The addition of oxygen yielded this new glow. Con- 
siderable further study of this glow is planned because of 
its probable importance for the understanding of excitation 
processes in the upper atmosphere. 


A2. Spectrophotometric Study of Reflectance Color 
Changes in Foods. E. J. Eastmonp, Western Regional Re- 
search Laboratory.—The influence of reflectance color on 
acceptability of foods has led to an investigation of the 
applicability of spectrophotometric methods to objective 
studies of color changes caused by processing conditions. 


Measurements have been made on a Hardy recording spec- 
trophotometer and on a Hunter multipurpose reflectometer. 
Various methods of sample preparation have been investi- 
gated and the precision of results compared. Measurements 
can be made on samples in the original form, but in many 
cases homogenization is necessary for precision. In certain 
dehydrated foods chromaticity has been found to be es- 
sentially constant over fairly wide ranges of particle size 
and moisture content even though the brightness may vary 
considerably. In other cases careful standardization of 
particle size would be necessary to measure small color 
differences. Attention has been given to tristimulus specifi- 
cation and the precision of determination of trichromatic 
coefficients. These coefficients are useful as indices of re- 
flectance color changes with time, temperature, and storage 
conditions. Luminous reflectance has also been found useful 
as a color variable. Changes in color are being correlated 
with quality and acceptability of food products. 


A3. Possibility of Altering the Decay Rate of a Radio- 
active Substance. Secré, University of California, 
Berkeley.—The radioactive decay constant of a substance 
decaying by orbital electron capture is proportional to 
|¥(0)|? of the electrons. In the case of a light element like 
Be’ it may be possible to alter this quantity by an appre- 
ciable amount by putting the Be in different chemical 
compounds. We would then have a slight change of the 
radioactive half-life of the Be in different compounds. The 
magnitude of the effect may be in the neighborhood of one 
percent, but it is practically impossible to give a quantita- 
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tive estimate because the total change of ¥(0) is affected 
by certain factors such as the density of the crystal, nature 
of the chemical bond, etc. They are both positive and nega- 
tive, and have comparable magnitudes. To obtain a reliable 
estimate of the effect we require a more detailed knowledge 
of the wave functions for various compounds than is at 
present available. Experiments are in progress to detect 
the effect by comparing the half-life of Be’ in Be metal with 
that in BeO or BeF 2. 


A4. Angular Distribution of e-Particles from Li’+ . R. 
F. Curisty AND S. Ruin, California Institute of Tech- 
nology.—One of us has measured the angular distribution 
of a-particles resulting from the bombardment of a thin 
target of Li by protons of from 0.4 to 1.4 Mev. The same 
experiments gave some information on the variation of the 
total cross section for the reaction in the same energy 
range. The angular distribution is well represented by 
1+A(E)cos*@. A increases with energy up to 0.9 Mev where 
it has a maximum value of 2.1 and then decreases slightly 
to a value of 1.8 at 1.4 Mev. An attempt has been made 
to understand in terms of the dispersion theory the angular 
distribution in conjunction with the magnitude and energy 
dependence of the total cross section. Although the cross 
section above 0.4 Mev is not well known, it seems to show 
no pronounced indication of resonance from 0 to 1.4 Mev. 
A reasonable fit of the data is obtained with two resonances 
with !'+1 Mev at.0 and 1.8 Mev roughly. The latter is 
taken to have J=2 but the former may be J=0 or J=2. 
The proton width agrees well with that derived from the 
y-ray resonance at 0.44 Mev only if the 0 energy resonance 
has also J =2. 


AS. Angular Distribution of Long-Range Alpha-Particles 
from Excited Ne**. SyLvAN RuBIN, California Institute of 
Technology.—Some preliminary measurements have been 
made of the angular distribution of the 5.9-cm alpha- 
particles from the reaction: F*(p, a)O"*. Measurements 
were made with a thin BeF; target at proton energies of 
1.14, 1.20, 1.33 Mev, corresponding to the two highest 
resonance peaks and the minimum between them on the 
excitation curve.! The angular distribution was measured 
from 25° to 155° to the proton beam. The distributions 
were markedly asymmetric, with much stronger yield back- 
ward than forward. Previous measurements near 400 kev? 


showed strong asymmetry also, but with greater yield: 


forward. Fitting of trial analytical angular distribution 
functions showed all powers of cos@, at least up to cos‘, 
are needed to fit the data. The cos*#@ and cos@ terms indicate 
an appreciable effect from d-wave incident protons. A large 
cos? term at 1.14 Mev, and large cos@ terms at 1.20 and 
1.33 Mev, indicate strong interference among adjacent 
long range alpha-emitting states of Ne. 
A. Pontes, and C. C. Pigs. Rev. 50, 258 
IW. B. McLean, A. Ellet, and J. A. Jacobs, Phys. Rev. 58, 500 (1940). 


A6. Device for Introducing Short-lived Radioactive 
Samples into a Cloud Chamber. T. Lauritsen, W. A. 
FowLer, anp C. C. LauriTsEN, California Institute of 


Technology.—In connection with the study of the radio- 
active products of Li* (0.88 sec.) a mechanism has been 
developed for automatically moving the activated target 
from the bombarding tube (at high vacuum) into the cloud 
chamber on a repeated cycle. The target is mounted at the 
end of a 3” stainless steel rod which, with an interlocking 
extension rod, passes through a series of locks connecting 
the bombardment position and the cloud chamber. Im- 
mediately after bombardment of the target, an air cylinder 
advances the rod, pushing the target into a socket in the 
end of the extension rod and passing the assembly succes- 
sively through a rough vacuum region, a region at cloud 
chamber pressure, and into the cloud chamber. Once the 
target is in the cloud chamber, a second air cylinder retracts 
the extension rod, exposing the target. The reverse motion 
consists in covering the target, passing through the locks 
and uncovering in the bombardment position. Seals be- 
tween the pressure locks are provided by Koroseal sleeves 
lubricated with high vacuum oil. A slight rotation of the 
rod facilitates passage through the seals. The elapsed time 
for a 10” motion of the target into the chamber is about 
1.5 sec. A somewhat longer interval is allowed on the return 
stroke to permit pumping out in the rough vacuum lock. 
This work was carried on under contract with the Office of 
Naval Research. 


A7. Auger Showers. M. M. Mixts anp R. F. Curisty, 
California Institute of Technology.—Under the current as- 
sumptions of proton primaries in cosmic rays, the Auger 
showers (previously explained by primary electrons) de- 
mand a secondary source of energetic electrons or photons. 
A fast decaying mesotron (r~10~* sec.) has already been 
proposed to account for the bulk of the soft component and 
8-decay lifetimes. The object of the present calculation is 
to determine whether this mesotron is also sufficient to 
account for Auger showers and bursts in unshielded ioniza- 
tion chambers. The coincident burst data of Lewis, which 
disagreed violently with the primary electron hypothesis," 
has been examined in a preliminary way. If this experiment 
is to be explained at all in terms of ionization due to elec- 
trons, it will probably require initiating electrons of energy 
>10" ev produced predominately near the top of the 
atmosphere and with several electrons having considerable 
angular spread associated in one event. In terms of decay 
this would mean r=10-" sec. An alternative possibility 
that the bursts are due to nuclear disintegrations is being 
examined. An experiment with three colinear chambers is 
suggested. The center would not necessarily fire simul- 
taneously with the extremes according to these possibilities. 

1L. Wolfenstein, Phys. Rev. 67, 238 (1945). 


A8. The Interpretation of the East-West Effect. Ropert 
A. MILuikan, California Institute of Technology.—T. H. 
Johnson's equatorial measurements of the East-West effect 
indicate that the East-West dissymmetry is less at very 
high altitudes than at sea level. Suggestions are herewith 
presented which would limit the unbalance between posi- 
tives and negatives to a region relatively near the earth in 
comparison with its radius and thus avoid the assumption 
of any electrical unbalance throughout space. Such an un- 
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balance has been generally recognized for fifty years as a 
valid objection to any theory which requires the earth to 
be bombarded by high speed particles of one sign only 
rather than of both signs, or else of no sign. If the validity 
of the principle is accepted all interpretation of the East- 
West effect must in some way or other conform with this 
principle of electrical balance throughout interstellar space. 
The suggestions necessary to produce such conformity are 
here discussed. 


A9. Cosmic Ray Directional Counter for the Deter- 
mination of the Azimuthal Variation of Primary Cosmic 
Radiation at Fixed Zenith Angles. ALFREDO BANos, JR. 
AND M. L. Perusguia, University of California, Los Angeles, 
and National University of Mexico, Mexico City.—At the 
1939 Symposium on Cosmic Rays held at the University 
of Chicago, M. S. Vallarta! suggested that the knowledge 
of the azimuthal variation of cosmic radiation could be 
used for the purpose of determining thesign and the energy 
spectrum of primary cosmic rays. In following this sugges- 
tion the writers designed, constructed, and operated, at the 
Institute of Physics of the National University of Mexico, 
an instrument consisting of four triple coincidence Geiger- 
Miiller telescopes mounted at zenith angles of 0°, 20°, 40°, 
and 60°, and so arranged that, by rotation around a vertical 
axis, the directional intensity or counting rate could be 
measured at 16 equally spaced azimuths. The rotation and 
photographic recording were made fully automatic. Several 
novel features of the apparatus are here discussed. Most 


important of all is the feature that, in contradistinction to 


the usual parallel connection of the three counters, as used 
in the standard Rossi* circuit and necessitating an amplifier 
tube for each Geiger-Miiller counter, a decided simplifica- 
tion in the circuits and a corresponding ease of operation 
were obtained by connecting the three counters in series 
feeding a single stage of amplification. Extensive statistical 
tests of the counting rate observed with the series con- 
nection were made by one of us* showing that the series 
connection was at least as reliable as the conventional 
parallel circuit. The results of two experiments, covering 
in excess of 100-days for each continuous run, have been 
reported by M. S. Vallarta.‘ 


1M. S. Vallarta, Rev. Mod. Phys. 11, 239 (1939). 
2B. Rossi, Nature 125, 636 (1930). 


* Alfredo Bafios, Jr., Rev. Mex. ge ing +, 19, 121-147 fp0et). 
S. Vallarta, M. L. Perusquia A de Bul Am. 


4M. 
Phys. Soc. 21, Cl (September, 1946). 


Al0. On the Meson Theory of Nuclear Forces. LEsLIc 
L. Fotpy, Department of Physics and Radiation Laboratory, 
University of California, Berkeley, California.—The failure 
of the Mgller-Rosenfeld-Schwinger mixture to yield quan- 
titative agreement with the observed properties of the 
deuteron has suggested the possibility that this agreement 
might be obtained with other mixtures in which tensor- 
force singularities are cancelled. Several such possibilities 
exist involving vector, pseudovector, and pseudoscalar 
mesons, some of which possess the interesting characteristic 
that they yield a finite tensor and central interaction in 
the limit where the meson masses approach equality and 
the square of the coupling constant increases indefinitely 


inversely as the difference in masses. In the . 
theories the resulting of the 


which is the same as would be obtained by differentiating 
the interaction with respect to the meson mass in the 
simple theories. The radial dependence of the central inter- 
action is attractive at short ranges but repulsive at long 
ranges and thus would presumably lead to preferential 
forward scattering of high energy neutrons scattered by 
protons. In spite of the admitted conceptual difficulties 
concerning the limiting process by which (1) was obtained 
and the field-theoretic significance of the result, the sim- 
plicity of the result was sufficiently provoking to induce an 
attempt to compare the consequences of (1) with experi- 
ment. Preliminary results appear to indicate that agree- 
ment with the properties of both the proton-neutron and 
proton-proton system cannot be obtained if one employs 
in (1) the mass of cosmic-ray mesons. 


All. Possible Use of Thermal Noise for Low Tempera- 
ture Thermometry. Epwarp GERJUOY AND A. THEODORE 
FORRESTER, University of Southern California.—The sug- 
gestion has recently been made that the Brownian motion 
of a quartz crystal be used to establish a low temperature 
thermodynamic scale.! The mean square noise voltage 
generated in any passive electrical network is 4kTR(f)df, 
where R(f) is the real part of the complex impedance. The 
noise voltage of a piezoelectric crystal, as well as that of 
ordinary resistors, must satisfy this relationship. The high 
impedance of a parallel resonant circuit at resonance sug- 
gests its use as a source of thermal noise. Subject to the 
condition that the output noise be representative of the 
low temperature to within 1 percent, a calculation was 
made which indicated that the lowest temperature measur- 
able in this manner is certainly greater than 0.5°K. This 
limitation arises from the fact that the input impedance of 
the first tube must be much larger than the impedance of 
the resonant circuit. As a result there is an upper limit to 
the impedance at resonance which does not permit the 
noise in the tube to be exceeded at very low temperatures. 
This consideration applies equally well to piezoelectric 
crystals as to ordinary circuits. 


1A, W. Lawson and E. A. Long, Phys. Rev. 70, 220 (1946). 


Al2. A Sensitive High Speed Radiation Thermocouple. 
Howarp Cary* anp K. P. GeorGe, National Technical 
Laboratories.—The case of a thermocouple receiving modu- 
lated radiation is analyzed to determine the conditions 
yielding the maximum signal to noise ratio. It is concluded 
that: (2) vacuum operation is best, although the advantage 
is small at high modulation frequencies; (b) the optimum 
electrical resistance decreases with increasing frequency; 
(c) the thermal capacity of the thermoelectric leads must 
be made negligible by making the lead length less than a 
limiting value which varies inversely as the square root of 
the modulation frequency. A vacuum thermocouple de- 
signed for optimum performance at,10 c.p.s. is described. 


ees 
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The receiver is of gold foil 2.0 mm by 0.3 mm by 0.2y thick, 
with the equivalent of 0.24 of gold. Taylor 

bismuth and bismuth-tin thermocouple leads 0.2 

mm long are used. To minimize thermal resistance of the 


receiver, four leads are soldered to the receiver at different 


ints and connected to form two parallel thermojunctions 
with a combined resistance of one ohm. In agreement with 
theory, steady radiation of approximately 2X10~' watt 
chopped at 10 c.p.s. and falling on the receiver produces an 
output signal equal to the r.m.s. noise of the thermocouple 
in a 1 c.p.s. band width. - 


* Now at Applied Physics Corporation, Pasadena, California. 


A13. A Proposed Experiment on the Ponderomotive 
Forces Exerted by Microwaves on Dielectric Media. Otto 


Havrern, Los Alamos Laboratory.—Some formulations! of 
relativistic electrodynamics lead to the result that an 
electromagnetic wave propagating itself in a dielectric 
medium exerts a force density given by the self-explanatory 
formula 

at 


The smallness of f has so far been considered to prevent 
experimental verification. The author proposes an experi- 
ment in which the force of a standing microwave in a \/4 
layer be measured preferably with a static magnetic field 
superimposed in the direction of the r. f. magnetic field. A 
(resonant) quartz-plate may be used as a piezoelectric indi- 
cator. This experiment if successful would permit for the 
first time to measure the actual time dependent of Poyn- 
ting’s vector. 


1Cf., e.g., W. Pauli, Relativitaets-Theorie. 


Bl. An Intense Positive Ion Source for Solids. A. E. 
Suaw, Argonne National Laboratory.—A positive ion source 
for solids has been developed. A tungsten crucible contain- 
ing the solid material and heated by electron bombardment 
produced intense beams of singly charged ions. Pt-Ir (15 
percent) yielded 120uA, for 100 watts of electron power, 
La:0; yielded 212uA, for 105 watts, and Sm;0; yielded 
40uA for 150 watts. A mass spectrographic analysis of the 
positive ions from a tungsten crucible containing La,O; and 
UN revealed the presence of Lat, LaOt+t, W+, WO*, U*+, 
UN*, UO* and possibly UO,*. A similar analysis of the 
ions from Pt-Ir (15 percent) revealed the presence of all 
the singly charged isotopes of W, Pt, and Ir. It is antici- 
pated that this source will find application in the mass 
deposition of isotopes. For example, in 10 minutes a de- 
posit of 0.50 microgram of the completely resolved isotopes 
of samarium has been obtained. 


B2. The Determination of Particle Size by the Low 
Angle Scattering of X-Rays. A. E. Sita, Shell Development 
Company.—The particle sizes of samples of pure y-Al:O; 
heated at various temperatures from 650°C to 1000°C were 
determined by the low angle scattering of x-rays by using 
a Philips Geiger counter x-ray spectrometer. The mean 
crystallite sizes of the same materials were also determined 


Taste I. 
Line Broadeni E.B.T 
Calculated Calculated Measured 
Specific Ss 
Surface, Surface, Surface, 
Heat Treatment Meanas m*/g MeanD m?*/g 
650°C— 6 hr. 86A 194 97A 198 172 
900°C—36 hr. 143 156 93.5 
1000°C—12 hr. 129 129 142 132 78 


by measurement of line broadening. The particle size and 
crystallite size, as measured by the above methods respec- 
tively, agree within the limit of experimental error. The 
surface calculated from the particle size distribution and 
mean crystallite size is in reasonable agreement with the 
surface measured by the Brunauer, Emmet, Teller method. 
The results are summarized in Table I. 


B3. Temperature Coefficient of Electrical Resistivity for 
Crystalline Selenium Containing Various Percentages of 
Bromine. WayNE E. BLACKBURN, Westinghouse Research 
Laboratories and University of Pittsburgh.—The tempera- 
ture coefficient of electrical resistance of brominated se- 
lenium over the temperature range, —30°C to 100°C, was 
found to be negative and variable. The electrical resistivi- 
ties of the selenium samples varied by factors of three to 
twelve for the temperature range mentioned depending 
upon the bromine content. A reproducible hysteresis loop 
was obtained on the plot of resistivity versus temperature. 
A minimum in the plot of electrical resistivity versus 
bromine content was found with 0.007 percent bromine in 
selenium. Activation energies for freeing electrons obtained 
from the slopes of conductivity versus temperature plots 
indicate a dependence upon bromine content in the 
selenium. 


B4. Viscous and Thermal Absorption of Sound as Re- 
laxation Phenomena. Cart Ecxart, Marine Physical 
Laboratory.*—In discussing the absorption of sound, refer- 
ence is often made to relaxation phenomena, as dis- 
tinguished from the “‘classical”’ viscous and thermal dissi- 


pative processes. While this may be convenient, it is not a 


basic distinction. By using known formulae, it is shown that 
viscous dissipation of sound can be considered as a relaxa- 
tion phenomenon, with relaxation time r,=4(»+-»’)/3pc*, 
where »,»’ = ordinary and dilational coefficients of viscosity, 
p=density of fluid, and c= velocity of sound. Similarly, the 
relaxation time for thermal dissipation is r,=«/pc*S, where 
«= thermal conductivity, and S=specific heat at constant 
volume. In the case of an ideal monatomic gas, the ratio 
s=r,/r, has the value 1.05 and both r, and r, are of the 
order of magnitude of one mean free time. For actual gases, 
s ranges from 0.57 for argon to 7.2 for air. In the case of 
liquids, a normal value of s appears to be 10, but it ranges 
from 0.032 for mercury to 9700 for glycerol. For liquids, 
there seems to be no simple picture of the two relaxation 
times, and there is no adequate theory of the ratio s for 
either gases as liquids. 
contract with the Bureau of Shipe and Office of Naval Research, Navy 
ureau a’ Navy 
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BS. Ionization Currents in Cylindrical Electrode Sys- 
tems. G. W. Jounson, University of California, Berkeley.— 
Ionization currents were measured in pure hydrogen be- 
tween concentric cylinders for pressures from 0.01 mm to 
740 mm. Above 1 mm 1/ log i/is is a function of X/p. At 
lower pressures it depends also on the pressure. Auxiliary 
anodes permitted estimates of the location of ionization in 
the field which confirmed the spatial distribution predicted 
by Morton.' Electrons released in the low and collected in 
the high field region resulted in about 10 percent less ion- 
ization than in the reverse case. In both cases the ionization 
‘was greater than given by Townsend's equation, as noted 
also at high pressures by Fisher and Weissler* for confocal 
paraboloids. The peaks of Morton’s i/i vs. p curves are 


identified with the pressure corresponding to the minimum . 


sparking potential. Application of the back diffusion equa- 
tion of Rice* in the form i= AV/(V+B/p), where A and B 
are constants, permitted the evaluation of i» at higher 
pressures than hitherto has been possible. The present work 
makes possible the delineation of tables and curves to 
predict the electron multiplication in highly divergent fields. 
1 P. L. Morton, Rev. 70, 358 
- 66, 95 (1944). 


?L. H, Fisher and G. L. Weissler, 
+C. W. Rice, Phys. Rev. 70, 228 (1946). 


B6. Negative Corona in Freon-Air Mixtures at N.T.P. 
G. L. WEIssSLER, University of Southern California.—In a 
point-to-plane corona gap of 3.1 cm, the onset potentials 
of the intermittent corona region were observed to be ap- 
proximately constant, 6-7000 volts (X/p~100) for all 
freon air mixtures, whereas in pure freon this value changed 
abruptly to 12,650 volts (X/p~200). It is suggested that 
this is mainly due to a marked decrease in the first Town- 
send coefficient alpha, similar to the findings of Kruithoff 
and Penning! in A-Ne mixtures. Concentrations of freon of 
10~ to one percent were very similar to the known behavior 
of air,? mixtures from one to twenty percent showed in the 
intermittent corona region Trichel bursts, that indicated 
a very rapid decay in the amplitude of individual pulses. 
This is due to the formation of negative ions of freon de- 
composition products, such as chlorine and fluorine, which 
must be more stable at higher values of X/p (about 200) 
than those of oxygen (about 90), Also their probability of 
electron attachment is greater, since higher energy elec- 
trons can dissipate their excess in numerous vibrational 
states of freon, and in dissociation. Thus negative space 
charges must be stable much closer to the point and there- 
fore be more efficient in choking off the discharge (decaying 
pulse size). Still higher concentrations inhibit proper 
avalanche formation such that Trichel bursts are barely 
observable. Current vs. potential curves show rapidly de- 
creasing slopes with increasing freon content. 


1 Kruithoff and Penning, Physica 3, 515 (1936); 4, 430 (1937). 
2 Loeb, Kip, Hudson, and Bennett, Phos Rev. 60, 714 AY 


B7. Positive Corona in Freon Air Mixture at N.T.P. 
E. I. Monr, University of Southern California (introduced 
by G. L. Weissler).—In order to further our understanding 
of positive corona and thereby of spark breakdown, space 
charge phenomena and current vs. potential character- 
istics were studied in a 3.1-cm point-to-plane gap (point 


diameter = $ mm) in freon air mixtures. Streamer a 

and strength increased going from pure air te 1 percent 
freon. Here the visual appearance of the corona changed 
markedly indicating that the discharge maintained itself 
solely by the streamer mechanism.' This seems due to a 
large increase in efficiency of photo-ionization in the gas 
without any appreciable change in the first Townsend 
coefficientgx. As the freon content was increased to higher 
and higher values, the streamers became successively more 
degenerate until they could not be distinguished from burst 
pulses. This is compared to similar observations in air using 
needle points; here only incipient streamers can 
because of the rapid decrease of the field near the point 
resulting in too small values of /2,°« dx. In our case the 
same was true because of the increased stability of negative 
freon decomposition product ions at higher values of X/p 
than those of Oz. The onset potentials of all observed 
phenomena were shifted to higher values and the current 
vs. potential curves showed decreasing slopes with increds- 
ing freon content. The changes in visual corona character. 
istics can be correlated to the space charge phenomena. 
The action of freon in suppressing discharge is ascribed to 
a gas mechanism, in contrast to the action of alcohol in 
A-counters,? which occurs at the cathode. 


1L. B. Loeb and A. F. Kip, J. Phys. 10, 142 (1939). 
2C. G. and D. D. . Frank. Inst. 447 (1941). 


B8. A New Derivation of the Method of Characteristics 
for Axially Symmetrical Supersonic Flow. A. Vazsonv1, 
North American Aviation, Inc. (introduced by S. H. Browne). 
—Two-dimensional supersonic flow in nozzles and around 
airfoils and other bodies can be treated by the “method of 
characteristics.” Recently this method has been extended 
to axially symmetrical flows. The differential equations 
using characteristic coordinates can be derived from the 
theory of hyperbolic partial differential equations. The 
derivation is laborious and does not clarify the basic 
physical principles involved. The alternative procedure 
proposed by the writer does not require an extensive knowl- 
edge of the theory of partial differential equations. Ber- 
noulli’s equation, the continuity equation, and the kine- 
matical expression for the rotation are written down using 
intrinsic coordinates (coordinates along streamlines and 
normal to streamlines). The differential change in velocity 
can be expressed as a function of the differential distance, 
the local value of the rotation, Mach number, distance 
from the axis of the flow, differential change in flow direc- 


tion, and the tangential and normal derivatives of the | 


velocity. The last two quantities drop out when changes 
along a Mach line are considered, and the basic equation of 
the method of characteristics is obtained. It is believed that 
this new derivation gives a better physical understanding 
of the method. 


B9. Atomization of Liquids. J. M. Scumipt, California 
Institute of Technology.—Various theories of the formation 
of droplets in pressure injection and air stream atomization 
of liquids are discussed. A number of experimental methods 
are described for the measurement of mean droplet sizes 
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as well as droplet size distributions in atomized liquids. 
Some of the methods used are high speed photography, 
impingement of the droplets on smoked plates, and photo- 
metric measurements of light absorption by sprays. Droplet 
size distributions are given for various liquids atomized 
under various conditions. In sprays of hydrocarbon fuels, 
injected under a pressure of 50 Ib. per in.* through hollow 
cone injectors into air streams moving with velocities of 
20 ft. per sec., mean droplet diameters as small as 45 mi- 
crons have been measured. 


B10. Escape from the Earth by Multiple-Step Rockets. 
MAarTIN SUMMERFIELD, California Institute of Technology. 
—Recent achievements in rocket vehicle performance as 
exemplified by the “V-2"’ rocket and the “WAC Corporal” 
rocket have focused attention on the physical factors limit- 
ing the velocities attainable in this manner. This paper 
deals with the question of escape from the surface of the 
earth, with special reference to four missions: (1) an earth- 
satellite revolving around the earth just above the atmos- 
phere, (2) a “stationary” earth-satellite whose angular 
velocity equals that of the earth, (3) a vehicle that com- 
pletely escapes the earth, (4) a vehicle that escapes the 
solar system. General performance equations are presented 
that bring out the effects of the following important 
parameters: jet exhaust velocity, ratio of propellant mass 
to gross mass, time of burning, and a drag factor. It is 
shown that, if exhaust velocities are limited by the energies 
produced by chemical reactions, it is doubtful that a single 
rocket can be constructed that is capable of leaving the 
earth. By means of multiple-step rockets much greater 
velocities can be imparted to the vehicle. Several cases of 
multiple-step escape rockets are calculated, and their 
characteristics discussed. Special emphasis is placed on the 
importance of determining the optimum number of steps. 


Bil. Problems in the Application of Nuclear Energy to 
Rocket Propulsion. H. S. Seirert anp M. M. MILLs, 
California Institute of Technology.—Consideration is given 
to the basic physical limitations involved in the application 
of nuclear energy sources to the propulsion of rockets. 


Quantitative estimates are made of the performance 
parameter specific impulse (impulse transmitted to the 
rocket per unit mass of expelled propellant) for the cases 
where propulsion is accomplished by: (a) photons, (b) 
fission fragments, and (c) an inert working fluid. Because 
of the fundamental difference between chemical rockets, in 
which the propellant is used both as a source of energy and 
as the mass to be ejected in obtaining thrust, and a 
“nuclear” rocket of type (c) above, in which these two 
functions may be separated, there exists an optimum ratio 
of payload or nonexpendable mass to working fluid or 
expendable mass, which differs for the two types of rockets. 
This was first pointed out by Ackeret.* 

* J. Ackeret, Helv. Phys. Acta 19, 103 (1946). 


B12. Gamma-Radiation from Be’+H'. C. C. LAuRITSEN, 
T. Lauritsen, AND W. A. Fow.er, California Institute of 
Technology.—Further investigations have been made of the 
gamma-radiation from Be® bombarded by protons! par- 
ticularly in the region of the resonances in the thin target 
excitation curve near 0.975 and 1.06 Mev. The first reso- 
nance has been found to have a width of 100 Kev. The 
thick target yield is 2 10-* quanta per proton, so that 
#I',T,/T equals 20 volts. The second resonance width is 
less than 10 Kev, higher resolution in proton energies and 
thinner targets being necessary to complete this determina- 
tion. The thick target yield is 2 10~* quanta per proton so 
that wI",T',/T equals 2 volts. Measurements of the second- 
ary pairs from a thick target bombarded at 1.1 Mev indi- 
cate radiation of energy 7.20.2 Mev. Absorption measure- 
ments indicate no soft radiation. Thus the radiation from 
the strong wide resonance, at least, arises from a single 
transition from the compound nucleus B" to its ground 
state. The width in this case must be due to competition 
with the re-emission of s-wave protons and probably with 
the emission of deuterons and alpha-particles. The narrow 
resonance must arise from p-wave capture with alpha and 
deuteron emission relatively improbable or forbidden. This 
work was carried on under contract with the Office of Naval 
Research. 

1 W. J. Hushley, Phys. Rev. 67, 34 (1945). 
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